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SECTION  I 
INTRODUCTION 

The  calculation  of  flow  parameters  in  a  supersonic  airstream  is 
accomplished  by  application  of  the  well-known  conservation  equations  of 
fluid  mechanics.  These  equations  basically  state  that  for  a  steadily 
flowing  (time  independent)  system,  the  mass,  momentum,  and  energy  leaving 
a  control  volume  must  be  equal  to  the  mass,  momentum,  and  energy  entering 
the  system,  assuming  that  there  is  no  change  in  these  quantities  within 
the  control  volume.  The  application  of  the  conservation  equations  to  the 
various  phenomena  that  can  occur  in  a  supersonic  airstream  is  quite  straight¬ 
forward,  and  for  low  supersonic  Mach  numbers  is  well  documented.  For 
higher  supersonic  Mach  numbers  where  the  ratio  of  specific  heats  (y) 
cannot  be  assumed  constant,  the  documentation  is  not  quite  as  complete. 

The  purpose  of  this  report  is  to  present  in  a  concise  manner  the 
calculation  techniques  for  both  low  and  high  supersonic  Mach  number  air¬ 
flows.  The  flow  phenomena  addressed  include  normal  shock  waves,  oblique 
shock  waves,  and  isentropic  expansions.  While  all  the  techniques  included 
are  generally  applicable  only  to  two-dimensional  flow,  the  special  case 
of  axisymmetric  flow  over  a  cone  at  zero  angle  of  attack  is  also  addressed. 

In  addition,  oblique  shocks  resulting  from  flow  over  swept  wings  are 
considered. 

The  calculation  techniques  for  low  supersonic  Mach  numbers  are 
presented  in  Section  II.  These  techniques  employ  constant  y  equations 
and  are  quite  straightforward.  Separate  subsections  under  Section  II  cover 
normal  shock  waves,  two-dimensional  oblique  shock  waves,  isentropic 
expansions,  and  oblique  shock  waves  for  cones  at  zero  degrees  angle  of 
attack.  Included  in  each  subsection  is  a  discussion  of  a  computer 
subroutine  that  has  been  written  to  be  used  as  a  computation  aid. 

Section  II. 1  covers  two-dimensional  oblique  shock  waves.  Included 
in  Section  II. 1  is  the  derivation  of  a  closed  form  solution  for  deter¬ 
mining  shock  wave  angle,  given  the  upstream  Mach  number  and  the  deflec¬ 
tion  angle.  An  additional  feature  of  Section  II. 1  is  the  inclusion 
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of  charts  that  can  be  used  in  determining  properties  downstream  of  an 
oblique  shock.  These  charts  have  been  drawn  to  include  upstream  Mach 
numbers  as  high  as  Mach  4.  Section  II. 4  (on  shock  waves  emanating  from 
the  apex  of  a  zero  degree  angle  of  attack  cone)  also  includes  charts 
useful  for  determining  conditions  downstream  of  the  shock  wave.  Example 
cases  illustrating  the  use  of  the  charts  are  included. 

Calculation  techniques  for  high  Mach  number  supersonic  flows  are 
presented  in  Section  III.  As  in  Section  II,  separate  subsections  are 
included  for  each  of  the  various  flow  phenomena.  For  the  high  supersonic 
Mach  number  solutions,  precise  closed  form  solutions  are  not  available, 
and  the  calculation  techniques  are  all  iterative  solutions  employing  the 
conservation  equations  and  an  approximate  model  for  equilibrium  air 
devised  by  Hansen  and  Hodge  as  reported  in  NASA  TN  D-352.  A  versatile 
computer  subroutine  based  on  this  model  has  been  developed  and  is 
described  in  Section  III.  Each  of  the  computer  subroutines  written  for 
the  various  flow  phenomena  employs  this  subroutine  for  equilibrium  air 
calculations.  Because  of  the  number  of  variables  involved,  a  set  of 
charts  such  as  presented  in  Section  II  is  not  possible  for  calculation  of 
the  high  Mach  number  flow  parameters.  The  computer  subroutines  presented 
in  Section  III  are  the  most  expedient  means  for  obtaining  the  required 
information. 

Computer  programs  employing  the  subroutines  are  included  in  Section  IV 
for  several  hypothetical  cases.  Solutions  are  also  worked  out  for  each 
of  the  cases,  using  the  charts  described  in  Section  II,  in  order  to 
illustrate  the  magnitude  of  the  error  resulting  from  the  use  of  the 
charts  at  various  Mach  numbers.  All  of  the  examples  discussed  in  Section 
IV  are  actually  very  idealized  cases.  No  flow  process  is  truely  two- 
dimensional;  since  interactions  with  sideplates  or  edges  will  occur  in 
flow  over  any  surface  of  non-infinite  width.  Additionally,  the  effects 
of  the  viscosity  of  the  air  have  not  been  considered  in  the  calculation 
techniques.  Viscous  interactions  and  three-dimensional  effects  are 
beyond  the  scope  and  purpose  of  this  report.  Since  these  effects  that 
occur  in  any  real  situation,  are  not  included,  the  results  will  not 
necessarily  be  in  close  agreement  with  data  obtained  from  actual  test 
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hardware.  To  predict  data  as  obtained  from  test  hardware,  one  must 
employ  empirical  methods,  derived  from  data  obtained  at  similar  conditions, 
in  conjunction  with  the  calculation  techniques  reported  herein.  This 
also  is  beyond  the  scope  and  purpose  of  this  report.  Therefore,  no 
comparisons  with  test  data  are  included  in  the  report. 
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SECTION  II 

TECHNIQUES  FOR  LOW  MACH  NUMBER  (M<4 )  SUPERSONIC  AIR  FLOWS 

The  well  known  calculation  techniques  described  in  this  section  are 
included  primarily  for  completeness.  For  each  of  the  flow  phenomena 
considered  (normal  shock,  oblique  shock,  isentropic  expansion,  and  conical 
shock),  a  computer  subroutine  was  written  to  perform  the  required  cal¬ 
culations.  These  subroutines  are  available  for  use  in  a  computer  program 
that  can  easily  be  written  to  handle  any  two  dimensional  flow  situation. 
They  were  also  employed  to  calculate  the  data  plotted  on  the  charts 
associated  with  each  of  the  calculation  techniques. 

1 .  NORMAL  SHOCK  WAVE 


S77777777777777777777 


For  flow  calculations  involving  a  normal  shock,  the  Mach  number, 
static  pressure,  and  static  temperature  upstream  of  the  shock  wave  are 
generally  known,  and  it  is  desired  to  find  the  value  of  those  parameters 
downstream.  Also  of  interest  is  the  entropy  increase  or  total  pressure 
ratio  across  the  shock.  Equations  listed  in  NACA  Report  1135  can  be 
used  to  calculate  these  quantities.  The  derivation  of  the  equations  is 
trivial  and  will  not  be  repeated  here.  The  equations  to  be  used  are 


P2  2  y  M^  -  (y-1) 

P-|  Y  +  1 

T2  [2  Y  -  (Y-1)]  [(y-1)  +  2] 

Ti  (y+D2Mi 


(1) 


(2) 


M 


2 


/(y-D  M 


2  y  Mf  - 


+  2 
(y-D 


(3) 
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For  convenience,  the  normal  shock  subroutine  also  includes  the 
calculation  of  total  pressure  and  temperature  upstream  and  downstream  of 
the  shock  wave.  These  parameters  are  found  by  applying  the  following 
equations : 

-y  =  1  +  M2  (5) 


Using  Equations  1  through  6,the  desired  quantities  may  be  easily 
calculated  using  non-iterative  methods.  The  calling  statement  for  the 
computer  subroutine  is 

CALL  NORM  (GAM,  T1 ,  PI,  Ml,  T2,  P2,  M2,  PR,  TT,  PT1  ) 

The  first  four  items  in  the  sequence  are  inputs.  The  remaining 

items  are  outputs,  and  for  the  most  part  are  self  explanatory.  The 

units  of  the  output  quantities  will  be  the  same  as  the  units  of  the 

corresponding  input  quantities.  The  quantity  PR  is  the  total  pressure 

ratio  (P.  / P.  ).  TT  is  the  total  temperature,  and  PTI  is  the  upstream 
l2  L1 

total  pressure.  The  analyses  herein  are  based  on  the  assumption  that  the 
system  being  considered  is  adiabatic;  thus,  the  total  temperature 
upstream  and  downstream  of  the  shock  is  the  same.  The  entropy  change 
across  the  shock  is  related  to  the  total  pressure  ratio  by  the  equation 
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Subroutine  NORM  (Table  A-1)  was  used  to  calculate  the  data  presented 
in  Figure  B-l.  Figure  B-l  can  be  used  in  performing  hand  calculations 
for  systems  containing  a  normal  shock.  Data  are  plotted  in  Figure  B-l 
for  Mach  numbers  up  to  4.  The  method  described  in  Section  III  should 
be  used  for  Mach  numbers  greater  than  about  4  and  for  cases  in  which 
the  flow  upstream  of  the  shock  wave  has  an  elevated  temperature. 

2.  OBLIQUE  SHOCK  WAVE-PLANAR 


A  supersonic  2-D  airflow  encountering  a  positive  deflection  will 
experience  an  oblique  shock  wave.  Generally  the  conditions  upstream 
of  the  shock  wave  and  the  flow  deflection  angle  are  known,  and  it  is 
desired  to  determine  the  conditions  downstream  and  the  angle  of  the 
shock  wave.  The  derivation  of  the  oblique  shock  equations  is  also  well 
documented  and  will  not  be  repeated  here.  The  equations  listed  in  NACA 
1135  (Reference  2)  for  oblique  shock  waves  will  be  used  in  the  computer 
program.  The  equations  used  are 


2  y  M|  sin^  0  -  (y-1 ) 
y+1 


(8) 


[2yM2  sin2  0-(y-l)]  [(y-1)  M2  sin2  9+2] 

(y+1 )2  M2  sin2  0 

(y+1 )2  M^  sin2  0-4(M2  sin2  0-l)(yM2  sin2  0+1) 
[2yM2  sin2  0- (y-1 ) ][ (y-1 )  M2  sin2  9+2] 


(9) 


(10) 


% 

(y+1)  M2 

sin2  0 

y-1 

y+1 

= 

(y-l)  M2 

sin2  0+2 

2yM2  sin2  e-(y-l) 

1 

y-1 


(ID 
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Equations  8  through  11  can  be  easily  solved  if  the  shock  wave  angle 
is  known.  However,  as  stated  previously,  the  deflection  angle  is 
generally  the  known  parameter,  and  shock  wave  angle  is  one  of  the  para¬ 
meters  to  be  determined.  An  equation  relating  the  deflection  angle  and 
the  shock  wave  angle  is 


cot6 


tanG 


(y+1)  M? 

2(M^  si n^  6-1 ) 


-1 


(12) 


Examination  of  Equation  12  will  show  that  the  use  of  this  equation 
would  require  an  iterative  solution.  To  avoid  the  iterative  solution 
another  equation  from  NACA  1135  (Reference  2)  is  employed: 


where 


sin6  0 


psin  6  +  q  sin  0  +  r  =  0 


(13) 


P  =- 


Mi  +  2 
M? 


-  y  sin  6 


2Mj  +  1 


+  Izl 
2 

M1 


2 

sin  6 


(14) 


(15) 


m 


2„ 
cos  5 


(16) 


Solving  Equation  13  explicitly  for  0  is  quite  complicated,  so  the 


solution  will 

be  detailed  below. 

If  we  let: 

.  2 

y  =  sin  0 

(17) 

Then : 

3  2 

y  +  py  +qy  +  m  =  o 

(18) 

If  we  let: 

y  =  ( x  -  3 ) 

(19) 

Then : 

x^  +  b-|  x  +  b^  =  0 

(20) 
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Where:  b]  =  1/3  ( 3q-p2 )  (21) 

b2  =  1/27  (2p3-9pq  +  27m)  (22) 

The  solutions  to  Equation  20  are 


x 


k 


cos  (|  +  120°k  )  ,  k 


0,  1,  2 


(23) 


cos  4>  =  +^y-27b|/4b^  (24) 

In  Equation  24  the  upper  sign  is  used  if  bpis  positive,  and  the 
lower  sign  is  used  if  b^is  negative. 

The  strong  shock  solution  is  found  from  Equation  23  with  k=0.  The 
weak  shock  solution,  which  is  the  one  generally  of  interest,  is  obtained 
by  setting  k=2  in  Equation  24.  Using  this  value  of  x^  in  Equation  19  and 
then  solving  for  0  in  Equation  17  provides  the  shock  angle  of  interest. 

The  solution  obtained  with  k=l  is  of  no  physical  significance.  Once  the 
shock  angle  has  been  obtained,  it  is  a  simple  matter  to  determine  the 
other  parameters  of  interest  by  applying  Equations  8  through  10.  The 
calling  statement  for  the  oblique  shock  computer  subroutine  is 

CALL  CGSK  (DELT,  VM1 ,  T1 ,  PI,  GAM,  SQ,  DELTM,  THETA,  VM2 ,  T2 ,  P2) 

The  first  six  items  in  this  sequence  are  inputs.  DELT  is  the 
deflection  angle,  VM1  is  the  Mach  number  upstream  of  the  deflection,  and 
T1  and  PI  are  the  temperature  and  pressure  upstream  of  the  deflection. 

GAM  is  the  value  of  specific  heat  ratio  to  be  used  for  the  calculations. 

If  a  strong  shock  solution  is  desired,  SQ  is  set  equal  to  one.  If  SQ  is 
set  equal  to  zero,  the  weak  shock  solution  is  calculated.  DELTM  is  the 
maximum  deflection  angle  for  which  an  attached  oblique  shock  can  exist 
for  the  input  value  of  Ml.  If  DELTM  is  less  than  the  input  value  of 
DELT,  no  solution  exists,  and  all  of  the  output  parameters  will  be  set 
equal  to  zero.  THETA  is  the  shock  angle  calculated,  and  VM2,  T2,  and  P2 
are  the  Mach  number,  temperature  and  pressure  downstream  of  the  shock. 
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The  units  of  temperature  and  pressure  will  be  the  same  as  those  of  the 
inputs.  All  angles  are  input  and  output  in  degrees.  Neither  the 
total  pressure  ratio  nor  the  entropy  increase  are  directly  calculated  in 
this  subroutine.  However,  with  the  static  pressure  and  temperature  known 
on  both  sides  of  the  shock  wave,  the  entropy  increase  is  easily  determined 
by  use  of  the  following  equation: 

T  P 

AS  _  ,2  .2 

R  7^-T  £n  T-|  "  £n  P1  (25) 

Total  pressure  ratio  is  then  easily  determined  by  application  of  the 
inverse  of  Equation  7: 


Pt,  AS 

tH2  =  e"  R  (26) 

*1 

Subroutine  CGSK  (Table  A-2)  has  been  employed  in  calculating  the 
data  presented  in  Figures  B-2  through  B-8.  Figure  B-2  indicates  the 
maximum  deflection  angle  possible  for  a  given  Mach  number  preceding  the 
shock.  The  shock  angle,  and  Mach  number,  temperature,  and  pressure  can 
be  found  by  use  of  Figure  B-3  through  B-6.  The  entropy  increase  for 
small  deflection  angles  is  plotted  on  Figure  B-7,  and  Figure  B-8  has  the 
same  parameter  for  larger  deflection  angles. 

3.  OBLIQUE  SHOCK  WAVE  -  SWEPT  WING 


Swept  Wing  Viewed  from  Above 
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While  the  equations  and  computer  subroutine  discussed  in  Section  II. 2 
are  st-'ctly  applicable  only  to  two-dimensional  situations,  it  is  possible 
that  by  appropriate  manipulation  one  can  apply  them  to  some  problems  of 
flow  over  a  swept  wing.  Referring  to  the  diagram  on  the  preceding  page, 
two  cases  are  possible  with  swept  wings.  In  the  first  case  the  Mach  angle, 
PM,  is  greater  than  the  swept  angle,  A.  In  this  case  the  Mach  number  nor¬ 
mal  to  the  leading  edge  of  the  wing  will  be  subsonic,  and  the  oblique  shock 
equations  cannot  be  applied.  In  the  second  case,  is  less  than  A,  M 

°°n 

is  supersonic,  and  the  oblique  shock  equations  can  be  used  if  equivalent 
values  to  use  for  the  upstream  Mach  number  and  deflection  angle  can  be 
determined.  The  following  sketches  will  be  useful  for  determining  these 
equivalent  values. 

Wing  Viewed  Parallel  to  Flight  Direction 
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The  equivalent  upstream  Mach  number  (M  ),  can  be  found  from 

e 

Equation  27  which  is  derived  as  follows: 


Similarly,  the  equivalent  angle  of  attack,  a  ,  and  wing  wedge 
angle,  a  ,  can  be  found  as  follows: 

-  *  -1  (  M°°Sina  \ 

ae  ~  an  l  M_cos  a  sinA  / 

00 

<»> 

And: 


tan  cr 

e 


t 

SL  sin  A 


tan  a  =  t /SL 


The  equivalent  deflection  angle  is 

5e  =  ae  +  ae  (3°) 

Charts  for  determining  Mro  /M^,  a g  and  ag  are  included  as  Figures  B-9  and 

and  B-10.  Also,  a  computer  subroutine  was  written  to  determine  the 
equivalent  Mach  number  and  deflection  angle  for  input  to  subroutine  CGSK. 
The  calling  statement  for  the  subroutine  is 

CALL  SWEP  (LAM,  M,  ALPHA,  SIG,  MEQ,  ALPEQ,  SIGEQ) 
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The  first  four  quantities  are  inputs  to  the  subroutine,  and  the 
last  three  are  the  desired  outputs.  Note  that  the  ratio  of  specific 
heats  does  not  enter  into  the  calculations  of  this  section.  Therefore, 
the  procedures  developed  in  this  section  and  the  computer  subroutine  are 
applicable  to  either  low  or  high  Mach  number  supersonic  airflows. 

See  Table  A-3  for  the  listing  of  subroutine  SWEP. 

4.  ISENTROPIC  EXPANSION 


A  two  dimensional  airflow  encountering  a  negative  deflection  or 
corner  will,  in  an  inviscid  system,  experience  an  isentropic  expansion. 
The  expansion  is  governed  by  the  well-known  Prandtl -Meyer  Equation 

v  =  tan_iy^|  (M2-!)  -  tan-1  ^  (31) 

The  Prandtl -Meyer  angle,  v,  is  a  property  of  the  airflow  and  is  a 
function  only  of  the  Mach  number  (for  a  const  y  system).  The  relation¬ 
ship  between  Prandtl -Meyer  angle  upstream  and  downstream  of  a  corner  and 
the  deflection  angle  is 

v2  =  v1  -  (  -6)  (32) 

Therefore,  if  the  Mach  number  upstream  of  a  corner  is  known,  the 
corresponding  value  of  v  can  be  found.  Adding  to  this  quantity  the 
deflection  angle  results  in  the  value  of  v  downstream  of  the  deflection. 
From  this  the  corresponding  value  of  Mach  number  downstream  of  the 
deflection  can  be  obtained  through  an  iterative  solution  of  Equation  31. 
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Since  the  expansion  is  isentropic  and  adiabatic,  the  static  temperature 
and  pressure  can  be  found  from  the  isentropic  flow  equations 


T 

T 


2 

1 


1  + 


1  + 


(33) 


Since  the  flow  is  isentropic, 


And 


As 

R 


=  0 


1  .0 


(34) 


An  additional  parameter  generally  used  for  calculations  involving 
isentropic  turning  is  the  Mach  angle,  y.  The  Mach  angle  is  simply  a 
function  of  the  Mach  number  and  is  defined  by  the  equation 


p  =  sin'1  1  (35) 

The  computer  subroutine  for  isentropic  turns  includes  the  calculation 
of  Mach  angle  upstream  and  downstream  of  the  turn.  The  calling  statement 
for  this  subroutine  is 

CALL  EXN  (DELT,  VM1 ,  T1 ,  PI,  GAM,  VMU1 ,  VMUi,  VM2,  T2,  P2) 

The  input  quantities  for  this  subroutine  are:  DELT,  the  deflection 
angle  in  degrees;  VM1 ,  the  upstream  Mach  number;  T1 ,  the  upstream  static 
pressure  in  °R;  PI,  the  upstream  static  pressure  in  atmospheres;  and  GAM, 
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the  ratio  of  specific  heats.  For  the  subroutine  to  function  correctly, 
the  value  of  DELT  must  be  negative.  The  parameters  determined  by  use  of 
subroutine  EXN  are:  VMUl  ,  the  Mach  angle  upstream  of  the  deflection  in 
degrees;  VMU2,  the  Mach  angle  downstream  of  the  deflection,  in  degrees; 
VM2,  the  Mach  number  downstream;  and  T2  and  P2,  the  static  temperature 
and  pressure  downstream.  The  listing  of  subroutine  EXN  is  included  in 
Table  A-4.  A  subroutine  using  Newton's  Method  is  employed  to  facilitate 
the  iteration  procedure.  A  listing  of  this  program  (subroutine  NEWT) 
is  included  as  Table  A-5. 

Figures  B-ll  and  B- 1 2  are  included  to  aid  in  hand  calculations  of 
isentropic  turning  situations.  On  Figure  B-ll  the  Prandtl -Meyer  angle 
is  plotted  as  a  function  of  Mach  number.  On  Figure  B-12  the  parameters 
T/Tt  and  P/P^  are  plotted  as  functions  of  Mach  number.  Since  T  and  P^ 
are  constant  for  this  situation,  the  static  temperature  and  pressure 
downstream  of  the  deflection  can  be  found  once  the  downstream  Mach  number 
has  been  determined. 

5.  OBLIQUE  SHOCK  WAVE  -  AXISYMMETRIC 


The  flow  over  a  cone  at  zero  degrees  angle  of  attack  is  a  special 
case  of  two-dimensional  flow.  The  two  dimensions  used  to  describe 
the  position  of  a  point  in  space  are,  in  this  case,  the  conical  coordinate 
angle  of  the  point  and  the  distance  to  the  point  from  the  cone  tip.  The 
solution  for  this  special  case  was  developed  by  Taylor-Maccol 1 .  The 
Taylor-Maccol 1  Theory  is  described  in  numerous  references,  but  will  be 
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repeated  here  for  completeness.  The  following  sketch  illustrates  the 
parameters  involved  with  an  element  of  fluid  downstream  of  the  shock 
wave : 


For  a  steadily  flowing  system,  the  Continuity  Equation  for  the 
elament  of  fluid  is 


pu  [rde  (2nr  sin6  )]  +  t.v  [dr  (2irr  sin6  )] 


=  [pu  +  9  dr] [ ( r+dr )  d6r (2tt) ( r+dr)  sinSr] 


+  [pv  +  der][dr(2:tr)  sin  (0„  +  dttj] 


(36) 


Multiplying  this  out  and  dividing  by  2mdrd9,  we  have,  after 
rearrangement , 


pvr 


sin0p  sin(Qr+d0r) 


dQ 


dQ 


=  2pur  sin0  +  pudr  sin9  +  r2  sin  9 

r  r  9r  r 


+  r  dr  sine^  +  (dr)2  sine^,  +  -3-^Pv-^-  r  sin  (0r+d0r)  (37) 


pu  dr  sin  0r  ■*  0 

r  dr  sin  0  -*•  0 
ar  r 

(dr)2  sin  0  -  0 
3r  '  r 


As  dr  approaches  0: 
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And  as  d-j^  approaches  0: 


sin  vr-  sin  (ur  +  dGr ) 


sin  6  -  sin  co 
r  r 


co^-d'*'^  -  cos  sjx'cf^', 


d9 


sin  6  -  sin  G  -  cos  9  d  9 
_ r _ r _ r_ _ r 

dS 


=  -COS  0 


sin  (O  +  d6  )  sin  0 
r  r  r 


Equation  37  then  becomes 
2pur  sir,  9_  +  sin  0_  +  - 


r  3r 


r  .  aQ  r  sin  9r  +  pvr  cos  9r  =  0  (33) 


The  flow  in  the  region  between  the  shockwave  and  the  surface  is 
assumed  to  be  irrotational .  This  is  illustrated  by 


To  keep  the  crosshatched  element  from  rotating,  the  following 
equality  must  be  satisfied: 

v(rde  )  +  (u  +  d6  )  dr  =  u  dr  +  (r+dr)  d0„  (v  +  dr) 
r  r  r  o  r 

r 


Multiplying  this  out  and  dividing  by  drdB^,  we  have 


UL  =  r  iv  +  v  +  lv  dr 

30..  3r  3r 


(39) 


As  ar  approaches  zero: 


If  dr"° 
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Equation  39  then  becomes 


3u 

39. 


v  - 


The  condition  that  all  properties  are  constant  along 


that 


3r 


=  0 


Using  this.  Equation  38  becomes 


2pur  sin  9  +  r  sin  6r  +  pvr  cos  9  =  C 

r  r 


And,  Equation  39  becomes 


v  = 


du 

d6 


After  expansion.  Equation  41  becomes 


2  pur  sin  0..  +  p  -  r  sin  6„  +  v  — -  r  sin  0„  + 


de 


Dividing  by  prsin0  ,  we  have 


2u  +  W  +  p  He  +  v  cot  er  =  0 

r  r 


Or 


2utWr  +  v(5  %  +  cot  er)  =  0 


Using  Equation  42,  Equation  43  can  be  written 


d  u  du 

de2  d9r 


(  cot  9r  +  p  He  )  +  2u  =  0 


(40) 

rays  specifies 

(41) 

(42) 

vr  cos  0=0 
r 

(43) 

(44) 
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To  make  integration  of  Equation  44  possible,  it  is  necessary  to 
replace  the  expression  involving  c  by  one  involving  u  and  fj  only.  If 
we  specify  that  the  flow  between  the  shock  wave  and  the  surface  is 
isoenergetic,  we  can  use  the  energy  equation  to  accomplish  this.  The 
form  of  the  energy  equation  to  be  used  is 


(u* 


2.  1 
V  )  =  1 


2 

c 


(45) 


Where  c  is  the  maximum  velocity  attainable  by  expanding  the  flow 
to  absolute  zero  temperature.  This  can  be  expressed  by 


c 


2 

.  U  SURF  \ 
2gJ  ) 


2 

c 


2 gj  h 


SURF 


+  u 


SURF 


(46) 


Since  we  have  specified  that  the  flow  is  inviscid,  the  flow  between 
the  shock  wave  and  the  surface  will  also  be  isentropic.  Therefore 


And 


P 

P 


t 

Y 

t 


>  =  P-  (P^1) 

P  _Y 


(P^1) 


Using  this  in  Equation  45,  we  have 


/  y-1 ,  .  1  /  2  ,  2 .  12 

(p  )  +  y  (u  +  v  )  =  c 


Now,  differentiating  this  with  respect  to  9  ,  we  have 


(47) 
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Or,  after  rearrangement. 


Substituting  Equation  48  in  Equation  47,  we  have 


-u 


Which  can  be  rearranged  to 


1  dp 

P  dQr 


+ 


) 


Izl 

2 


(c‘ 


-u 


2 


(48) 


Now,  this  can  be  used  in  Equation  44  to  eliminate  the  dependence 
on  density  variation.  The  resulting  equation  is 


(49) 


The  derivative  with  respect  to  0r  of  Equation  42  provides 

dv  d^u 


Using  this,  and  Equation  42  in  Equation  49,  we  have 
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Which  becomes,  after  collecting  terms 


^  - y2'  1i1(c2-u2)l=  (y-1  )(c2-u2)u  +  -^(c2-u2)cotv  v-vuv2-  Yy-cot  rv' 

fit!  1—  J 


This  can  be  nondimensional i zed  with  c  to  yield 


d2(“)  (l-  ^)c  +  *T-  O'  ^)COt9r  (c)  -y(c)(c)  '  i1  cot  r  (c) 


1+1/ vV  -  XnL 


0-4) 


Equation  50  can  be  solved  by  numerical  integration  in  an  inverse 
manner.  The  integration  proceeds  from  the  cone  surface  to  the  shock 
wave,  and  provides  properties  of  the  airflow  on  rays  emanating  from  the 
cone  apex.  The  technique  requires  that  the  cone  half  angle  and  the 
velocity  at  the  cone  surface  are  known.  The  calculated  parameters  are 
then  the  upstream  Mach  number,  and  the  shock  wave  angle. 


At  the  cone  surface 


Using  numerical  integration  the  first  derivative  on  any  ray  can  be 
expressed  approximately  as 
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Now,  if  this  value  of  the  first  derivative  is  averaged  with  the 
preceding  value,  and  added  to  the  preceding  value  of  u/c,  we  can  find 
the  new  value  of  u/c. 


With  Equation  52,  all  that  is  required  to  determine  the  parameters 
on  a  particular  ray  is  data  on  the  preceding  ray.  Therefore,  it  is 
possible  to  proceed  from  the  ray  on  the  cone  surface  (where  the  data 
is  known)  to  the  ray  on  the  shock  wave. 

An  example  of  the  numerical  integration  procedure  will  illustrate  its 
features.  For  the  example,  the  following  data  will  be  set: 

6  =  10°,  =  .8,  y  =  1-4,  A0r  =  1°  =  .01745  rad. 


At  the  surface  the  first  derivative  is  zero  (from  Equation  42). 
The  second  derivative  is  found  from  Equation  50: 


Then  proceding  to  the  next  ray  (at  9=11°)  and,  using  Equation  51, 


=  0  +  .01 745 ( - 1  .6)  =  -.027925 
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And,  using  Equation  52, 

(  =  .8  +  .01745(0)  +  .t-01?45)-2-  (-1  .6)  =  .799756 


Again,  using  Equation  50, 


.2  u\  ? 

cj  =  .4( . 36039) ( . 799756 )  +  .2(. 360 39 ) (cot  11 °) (- .027925)-! ,4( .799756) (-.027925) 

d02  /  1  .2(-.027925)2  -.2(. 36039) 

L  '  0  q 

rl  _  .2 (cot  11°)(-. 027925) 

1 ,2(-.027925)2  -.2(. 36039) 


-1 .46304 


This  procedure  is  repeated  for  successive  rays  until  the  assumed 
shock  wave  location  is  passed.  For  each  ray  it  is  convenient  to  deter¬ 
mine  the  resultant  flow  angle,  velocity  and  Mach  number.  These  can  be 
found  by  use  of  Equations  53  through  55: 

4>  =  tan  ^  +  0^  (53) 


(54) 

(55) 


A  summary  of  the  calculated  data  for  rays  out  to  22°  is  presented 
in  Table  1  for  the  example  case. 
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TABLE  1 

SUMMARY  OF  CALCULATIONS  FOR  CONE  FLOW 
NUMERICAL  INTEGRATION  EXAMPLE 


V 

u/c 

v/c 

w/c 

V 

M 

10° 

.8 

0 

.8 

10° 

2.981424 

n° 

.799756 

-.027925 

.800243 

9.000° 

2.983942 

12° 

.799046 

-.053460 

.800832 

8.172° 

2.990061 

13° 

.797904 

-.0773589 

.801645 

7.462° 

2.998547 

14° 

.796356 

-.100101 

.802623 

6.836° 

3.008818 

15° 

.794418 

-.122021 

.803734 

6.268° 

3.020573 

16° 

. 792102 

-.143379 

.804974 

5.740° 

3.033787 

17° 

. 789416 

-.164404 

.806354 

5.236° 

3.048629 

18° 

. 786364 

-.  185340 

.807910 

4.738° 

3.065546 

19° 

.782944 

-.206522 

.809724 

4.220° 

3.085486 

20° 

.779147 

-.228584 

.811986 

3.650° 

3.110717 

21° 

.774941 

-.253496 

.815349 

2.886° 

3.149000 

22° 

.770382 

-.268989 

.815992 

2.753° 

3.156426 

Despite  the  fact  that  rays  out  to  ©r  =  22°  have  been  calculated,  we 
as  yet  have  no  indication  of  the  location  of  the  shock  wave.  Data 
calculated  for  rays  upstream  of  the  shock  wave  is  not  valid.  To  find 
the  shock  wave  angle,  we  use  the  fact  that  the  flow  deflection  through  a 
conical  shock  wave  is  the  same  as  that  for  an  oblique  shock  wave  with 
the  same  upstream  conditions.  Also,  the  ratio  of  upstream  Mach  number 
to  Mach  number  just  downstream  of  the  shock  wave  will  be  the  same  for 
both  conical  and  oblique  shock  waves.  Therefore,  it  is  possible  to  use 
Figure  B-4,  which  was  plotted  for  oblique  shock  waves,  to  determine  the 
Mach  number  preceding  the  conical  shock  wave.  Then  Figure  B-3  can  be 
used  to  determine  the  shock  angle  for  the  conical  system. 
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To  demonstrate  this,  the  example  case  as  summarized  in  Table  1  will 
be  continued.  To  obtain  a  result  of  higher  precision  than  the  calculation 
increment  of  A0  ,  a  trial  and  error  solution  will  be  employed.  As  a 
first  guess,  we  shall  assume  the  shock  wave  angle  to  be  19°.  For  this  ray 
(from  Table  1)  the  flow  deflection  angle  is  4.220°.  Using  this,  and 
Figure  B-3,  we  find,  with  6/0  =  0.1936,  that  =  3.53.  Now,  using  this, 
the  deflection  angle,  and  Figure  B-4,  we  find  that  =  0.935  (3.53)  = 
3.300.  Checking  this  value  with  the  value  of  M  tabulated  in  Table  1, 
we  see  that  it  is  too  high.  Therefore,  the  next  guess  for  the  shock 
wave  angle  is  20°.  The  results  using  this  guess  are:  M-|  =  3.335, 

=  3.135.  Again,  this  value  of  M2  is  too  high,  so  the  next  guess  for 
the  shock  wave  angle  is  21°.  Using  this  guess,  the  results  are: 

Mi  =  3.077,  M2  =  2.923.  This  value  of  M2  is  too  low,  so  we  know  that  the 
shock  wave  angle  lies  somewhere  between  20°  and  21°.  To  obtain  a  more 
precise  result,  we  could  repeat  the  numerical  integration  between  6^  =  20° 
and  0r  =  21°  with  a  smaller  step  size.  An  alternate  method  is  to  simply 
interpolate  to  obtain  the  shock  wave  angle.  Using  this  method  we  find: 


-  ?n  +  1(3.135-3.1107) 

0  (3.135-3.1107)  +  (3.149-2.923) 


20.1° 


M2  -  3.1107  ♦  (3.1490-3.1107)  [(313'5_3^j^-^(130;)9_2  923)]  .  3.1144 
M2  -  3. 335-(3.335-3.077)  ^(3.1 35- 3 ) 1 1 07 ) ”+  ( 3?^49-2 ,923)J  *  3i31° 


A  computer  program  written  to  perform  the  calculations  illustrated 
in  the  preceding  example  is  included  as  Table  A-6.  The  input  data  to 
this  program  are:  DELT,  the  cone  half  angle  in  degrees;  USC,  the  surface 
velocity,  normalized;  GAM,  the  ratio  of  specific  heats;  and,  DTHETA,  the 
calculation  interval  in  degrees.  The  output  for  one  set  of  input  data 
is  on  a  single  sheet.  A  sample  output  is  included  as  Figure  1.  The 
input  for  this  sample  is  the  same  as  for  the  preceding  example: 

DELT  =  10.,  USC  =  0.8,  GAM  =  1.4,  and  DTHETA  =1. 
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Figure  1.  Sample  Output  for  Program  CONE 
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The  effect  of  the  calculation  interval  has  been  investigated  for  a 
case  with  a  cone  half  angle  of  20°,  a  normalized  surface  velocity  of  0.6, 
and  y = 1 . 4 .  The  results  of  this  investigation  are  shown  in  Figure  2, 
where  each  of  the  calculated  parameters  is  plotted  versus  2.6  .  Obviously, 
the  smaller  the  calculation  interval,  the  more  nearly  correct  the  numerical 
integration  results  will  be.  However,  there  is  a  practical  limit,  even 
with  a  high-speed  computer,  to  the  size  of  the  interval.  On  Figure  2  a 
point  has  been  placed  on  each  of  the  lines  indicating  the  value  of 
that  will  result  in  an  error  of  0.1S.  Thus  for  calculation  of  shock  angle, 
a  A9  =  1°  will  probably  be  sufficient  in  most  cases.  However,  if  the 
value  of  one  of  the  other  parameters  is  needed  with  great  accuracy,  a 
value  of  A9r=0.5°  or  less  may  be  needed. 

In  most  cases,  the  data  available  is  the  upstream  flow  conditions 
and  the  cone  angle.  The  calculation  technique  described  is  an  inverse 
technique,  and  therefore  must  be  iterated  in  order  t...  accommodate  the  usu¬ 
ally  available  data.  A  computer  subroutine  has  been  written  to  accomplish 
this.  The  calling  statement  for  the  conical  shock  computer  subroutine  is: 

CALL  CONSK  (VM1 ,  DELT,  GAM,  THE,  THETAS,  DSR,  VM2C ,  PHIC,  PRC,  TRC , 

VM2S,  PRSHK,  TRSHK,  PHIS,  VM2SRF,  PRSRF,  TRSRF) 

The  first  four  items  in  this  sequence  are  inputs.  VM1  is  the 
upstream  Mach  number,  DELT  is  the  cone  half  angle  in  degrees,  and  GAM  is 
the  ratio  of  specific  heats.  The  fourth  input  has  been  included  to 
provide  data  on  one  specific  ray  lying  between  the  shock  wave  and  the 
surface.  THE  is  the  ray  angle,  in  degrees,  for  which  data  is  required. 

The  output  information  is  as  follows:  THETAS  is  the  shock  wave  angle, 

DSR  is  the  non-dimensional i zed  entropy  across  the  shock  wave,  VM2C  is  the 
Mach  number  on  the  specified  ray,  PHIC  is  the  flow  angle  on  the  specified 
ray,  PRC  is  the  ratio  of  the  static  pressure  on  the  specified  ray  to  the 
upstream  static  pressure,  TRC  is  the  ratio  of  the  static  temperature  on 
the  specified  ray  to  the  upstream  static  pressure,  VM2S  is  the  Mach  number 
just  downstream  of  the  shock  wave,  PRSHK  is  the  ratio  of  static  pressure 
across  the  shock  wave,  TRSHK  is  the  ratio  of  static  temperature  across 
the  shock  wave,  PHIS  is  the  flow  angle  just  downstream  of  the  shock  wave, 
VM2SRF  is  the  Mach  number  on  the  cone  surface,  PRSRF  is  the  ratio  of  the 
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Figure  2.  Effect  of  Calculation  Increment  on  Accuracy 
of  Cone  Flow  Calculation  Procedure 
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static  pressure  on  the  surface  to  the  upstream  static  pressure,  and  TRSRF 
is  the  ratio  of  the  static  temperature  on  the  surface  to  the  upstream 
static  temperature.  A  listing  of  subroutine  CONSK  is  included  as  Table  A-7. 
Subroutine  NEWT  is  also  required  in  conjunction  with  subroutine  CONSK. 
Subroutine  CONSK  has  been  employed  in  calculating  the  data  presented  in 
Figures  B-13  through  B-17.  Figure  B-13  can  be  used  to  find  the  shock 
wave  angle  if  the  upstream  Mach  number  and  cone  half  angle  are  known. 

Figure  B-14  presents  the  Mach  number  ratio  across  the  shock  wave  and  ratio 
of  surface  Mach  number  to  upstream  Mach  number.  Figures  B-15  and  B-16 
can  be  used  to  find  the  ratio  of  upstream  static  temperature  or  pressure 
to  the  corresponding  parameter  just  downstream  of  the  shock  wave,  or  on 
the  cone  surface.  Figure  B-17  presents  the  non-dimensionalized  entropy 
increase  across  the  shock  wave.  The  data  plotted  in  these  charts  was 
calculated  with  a  calculation  interval  of  0.1  degree.  The  results,  then, 
are  quite  precise. 


To  illustrate  the  use  of  subroutine  CONSK  an  additional  program  has 
been  written.  The  listing  of  this  program,  which  is  entitled  program 
CONES,  is  included  as  Table  A-8.  The  input  data  to  this  program  are: 

VM1 ,  the  upstream  Mach  number;  DELT,  the  cone  half  angle  in  degrees;  GAM, 
the  ratio  of  specific  heats;  and  THE,  the  angle  in  degrees  of  a  ray  for 
which  data  is  desired.  A  sample  output  for  program  CONES  is  included  as 
Figure  3.  The  input  data  used  to  generate  Figure  3  is  as  follows: 

VM1 =3.0,  DELT=20 . ,  GAM=1.4,  THE=25. 
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Figure  3 
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SECTION  III 

TECHNIQUES  FOR  HIGH  MACH  NUMBER  (M>4)  SUPERSONIC  AIRFLOWS 

At  flight  Mach  numbers  greater  than  about  Mach  4,  the  temperature  of 
the  air  being  decelerated  by  a  propulsion  system  inlet  is  increased  to 
such  an  extent  that  the  equations  utilized  in  the  techniques  described 
in  Section  II  are  no  longer  applicable.  The  Section  II  equations  are 
generally  based  on  the  assumption  that  the  air  is  calorically  perfect, 
and  the  specific  heats  are  constant.  In  this  Section  the  conservation 
equations  will  be  employed  to  provide  "real  gas"  techniques  for  deter¬ 
mining  conditions  downstream  of  the  various  flow  pehnomena. 

The  real  gas  calculation  techniques  generally  require  iterative 
solutions  and  are  more  complicated  than  those  discussed  in  Section  II. 
Computer  solutions  are  required  for  all  but  a  few  sample  cases.  It 
is  not  possible  to  prepare  charts  for  use  as  computation  aids,  as  were 
presented  for  the  Section  II  techniques. 

Computer  subroutines  are  presented  for  each  of  the  flow  phenomena 
discussed  in  this  Section.  All  of  these  subroutines  employ  an  additional 
subroutine  that  provides  the  equilibrium  thermodynamic  properties  of 
air.  This  subroutine  is  based  on  the  approximate  model  for  equilibrium 
air  presented  by  Hansen  and  Hodge  in  Reference  1.  A  listing  of  subroutine 
THAIR  is  included  as  Table  A-9.  The  equations  used  in  this  subroutine 
will  not  be  repeated  here  since  they  are  taken  directly  from  Reference  1. 
The  use  of  the  subroutine,  however,  will  be  discussed.  The  calling 
statement  for  the  equilibrium  air  model  subroutine  is 

CALL  THAIR  (T,  P,  H,  S,  RHO,  AM,  CP,  CV,  GAM,  A,  K1 ,  K2 ) 

In  this  sequence: 

T  is  temperature  in  °R. 

P  is  pressure  in  atmospheres. 

H  is  enthalpy  in  Btu  per  pound. 

S  is  entropy  in  Btu  per  pound  °R. 

RHO  is  density  in  pounds  per  cubic  foot. 

AM  is  molecular  weight. 
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CP  is  the  specific  heat  at  constant  pressure  in  Btu  per  pound  -R. 

C V  is  the  specific  heat  at  constant  volume  in  Btu  per  pound  °R. 

GAM  is  the  local  isentropic  exponent,  defined  as  the  partial 
derivative  of  the  logarithm  of  pressure  with  respect  to  the  logarithm 
of  density,  at  constant  entropy.  In  equation  form,  this  is 

y  =  /  1  JP3-P-  \ 
rI  \ 3  log  o  /s 

(At  temperatures  less  than  about  1000°R,  the  isentropic  exponent  and  the 
ratio  of  specific  heats  have  very  nearly  equal  values). 

A  is  the  speed  of  sound  in  feet  per  second. 

K1  and  K2  are  operational  inputs. 

Subroutine  THAIR  has  a  number  of  options,  exercised  by  the  values 
assigned  to  K1  and  K2.  The  parameter  K1  may  be  assigned  a  value  from 
0  to  4.  The  effect  of  each  value  of  K1  can  be  shown  in  tabular  form: 

If  K1  =0  ■+•  Input  T  and  P;  Calculate  H,  S,  RHO,  AM 

If  K1  =1  -»•  Input  H  and  P;  Calculate  T,  S.  RHO,  AM 

If  Kl=2  -*■  Input  S  and  P;  Calculate  T,  H,  RHO,  AM 
If  K1  =3  -*■  Input  S  and  T;  Calculate  P,  H,  RHO,  AM 

If  K1  =4  -*  Input  S  and  H;  Calculate  T,  P,  RHO,  AM 

The  parameter  K2  may  be  assigned  a  value  of  zero,  one,  or  two.  The 
use  of  the  default  condition,  K2=0,  results  in  the  calculations  dependent 
on  the  value  of  K1 ,  as  listed  above,  and  values  are  not  determined  for 
CP,  C V,  GAM,  and  A.  If  K2=l ,  the  calculations  dependent  on  the  value  of 

K1  are  made,  and  in  addition  the  values  of  CP,  CV,  GAM,  and  A  are  deter¬ 
mined.  with  K2=2,  the  value  of  GAM  must  be  an  input  quantity,  and, 

depending  on  the  value  of  K1 ,  all  the  other  quantities  are  determined  on 

a  constant  gamma  basis.  This  option  is  particularly  useful  when  developing 
or  checking  out  a  new  computer  program  or  subroutine.  The  use  of  the 
K2=2  option  dramatically  reduces  the  computer  time  used.  The  K2 =1  option 
is  also  provided  to  reduce  computer  time.  The  calculation  of  CP,  CV,  GAM, 
and  A  takes  a  significant  amount  of  computer  time.  If  these  quantities  are 
not  required  for  a  particular  calculation,  the  computer  time  can  be  reduced 
by  setting  K2=0  and,  thereby,  eliminating  the  additional  computation. 
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More  precise  programs  for  determining  the  thermodynamic  properties 
of  air  are  available.  However,  the  use  of  these  programs,  particularly 
for  iterative  solutions,  can  be  quite  costly  in  computer  time.  Sub¬ 
routine  THAIR  is  a  good  compromise  between  precision  and  cost.  Comparison 
of  data  generated  by  subroutine  THAIR  and  data  in  standard  references, 
such  as  Reference  6  indicates  very  good  agreement  at  temperatures  up  to 
6000°R.  This  should  be  adequate  for  most  calculations  of  practical 
interest. 

1.  NORMAL  SHOCK  WAVE 

/  /  /  / 


For  any  flow  process  the  mass,  momentum,  and  energy  must  be  conserved. 
If  the  flow  is  steady  and  adiabatic  the  conservation  equations  for  the 
normal  shock  process  shown  above  can  be  expressed  as  follows: 


Energy  Conservation: 


(56) 


p2  v2  P]  V2 

Momentum  Conservation:  +  — - —  =  +  — - — 


(57) 


Mass  Conservation:  p2  V2  =  p-|  V-| 


(58) 


The  number  of  equations  can  be  increased  to  the  number  of  unknown 
quantities  by  applying  the  definition  of  Mach  number,  and  the  thermodynami c 
properties  of  air  as  follows: 


M 


2 


M 


(59) 

(60) 
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h1  =  ;•  (Tr  Pj) 

(61  ) 

a1  =  ;•  (T-,,  P^ 

(62) 

p  1  =  J  ( T 1  ,  p  1  ) 

(63) 

h  2  =  **  ^2 *  ^ 

(64) 

=  f  ( T 2 »  P 2 ^ 

(65) 

p2  =  (^2 5  ^2  ^ 

(66) 

Generally,  M-,  ,  T1  ,  and  P1  upstream  of  the  shock  wave  will  be  known, 
and  the  corresponding  properties  downstream  of  the  shock  wave  will  need 
to  be  calculated.  The  calculation  procedure  to  be  used  in  this  case 
is  as  follows: 

Step  1  -  Find  h-| ,  a-j ,  p-j  using  Equations  61,  62,  and  63 

Step  2  -  Find  V-|  using  Equation  60 

Step  3  -  Assume  a  value  for  p2 

Step  4  -  Find  V2  using  Equation  58 

Step  5  -  Find  P2  using  Equation  57 

Step  6  -  Find  h2  using  Equation  56 

Step  7  -  Check  the  assumed  value  of  p2  by  use  of  Equation  66 

Step  8  -  If  the  calculated  value  of  p2  (Step  7)  is  not  close 

enough  to  the  assumed  value  of  p2  (Step  3),  assume  a  new  value  of  p2  and 
repeat  Steps  4-7  until  the  assumed  and  calculated  values  are  close 
enough.  Then  proceed  to  Step  9. 

Step  9  -  Find  T2  and  a2  using  Equations  64  and  65 

Step  10  -  Find  M2  using  Equation  59 

Steps  1,  7,  and  9  will  require  a  Mollier  Diagram  for  equilibrium 
air  if  hand  calculations  are  performed  with  this  calculation  procedure. 
Such  a  diagram  is  contained  in  Reference  6. 

A  subroutine  for  calculating  on  a  real  gas  basis  the  properties 
downstream  of  a  normal  shock  has  been  written.  The  calling  statement 

for  this  subroutine  is: 

CALL  NRMSK  (VM1,  T1 ,  PI  A,  VM2 ,  T2,  P2A,  GAM) 
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The  input  items  for  this  subroutine  are  VM1  ,  T1 ,  P1A,  and  GAM.  VM1 
is  the  upstream  Mach  number,  T1  is  the  upstream  static  pressure  in  :R, 
and  P1A  is  the  upstream  static  pressure  in  atmospheres.  The  value  of 
GAM  is  zero  if  the  real  gas  solution  is  desired.  If  a  non-zero  value 
of  GAM  is  input,  the  calculations  will  be  done  on  a  perfect  gas  basis 
using  the  input  value  of  GAM  as  the  ratio  of  specific  heats.  This 
option  is  included  for  conservation  of  computer  time  while  developing 
new  programs  and  subroutines.  VM2,  T2,  and  P2A  are  the  calculated  para¬ 
meters,  and  have  the  same  units  as  the  corresponding  input  parameters. 

A  listing  of  subroutine  NRMSK  is  included  as  Table  A- 1 0 .  The  use  of 
subroutine  THAIR  is  required  with  subroutine  NRMSK. 

2.  OBLIQUE  SHOCK  WAVE-PLANAR 


The  problem  of  supersonic  flow  through  a  two-dimensional  deflection 
can  be  addressed  by  applying  the  conservation  equations  in  a  direction 
perpendicular  to  the  shock  wave.  A  vector  diagram  is  helpful  for 
illustrating  this: 


Then,  assuming  an  adiabatic  system,  the  conservation  equations 
perpendicular  to  the  shock  wave  are 


h 


2 


2 

U1 
2  g«J 


(67) 
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2  g  1  g 


P2U2 


From  the  vector  diagram: 

u.]  =  V1  sin  9 

The  upstream  Mach  number  is  related  to  the  velocity  by 


Where  the  speed  of  sound  is 

ai  Ri  Ti 

Combining  Equations  70  through  72: 

u?  -  M2  sin2  0(gY-j  R]  ) 

Equation  67  can  be 

/ 


Then,  substituting 


arranged  to 


with  Equation  73: 


M2  sin2  e(gv1  R-|  T] ) 
2gJ 


(68) 

(69) 

(70) 

(71) 

(72) 

(73) 


(74) 
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Although  this  is  a  real-gas  solution,  one  can  validly  use  relations 
involving  the  specific  heat  and  the  isentropic  exponent  as  long  as  the 
use  is  restricted  to  one  side  or  the  other  of  the  shock  wave.  The  use 
of  similar  relations  across  the  shocr.  wave  would  not  be  valid.  The 
upstream  temperature  and  enthalpy  are  therefore  related  by 


Where 


Equation  74  can  then  be  written 


Now,  rearranging  Equation  68: 


f  ?2  \  plul  '  P2U2 

h(pr’)  • 
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The  Perfect  3as  Law  is  valid  for  use  as  long  as  the  application  is 
for  conditions  on  only  one  side  of  the  shock  wave.  Therefore, 

P1  =  P1  R1T1 

Using  this  in  Equation  76: 


P 

P 


2 

1 


sin2  0(Yl) 


u, 

prp2  ~i 


u 


Using  Equation  69  (the  Conservation  of  Mass  equation): 


Which  becomes 


2  ?  ? 
pj-  -  1  =  Yt  si n  0 


Again  referring  to  the  vector  diagram  on  page  34: 


tan  (0-6) 


And 


1 

tan  0  =  — L 

y 


tan  (0-6)  =  ^2 
tan  0  u.j 


(77) 


(78) 
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Additional  equations  that  are  required  to  match  the  number  of 
unknowns  can  be  obtained  by  application  of  the  definition  of  Mach  number, 
the  vector  diagram  on  page  34,  and  the  thermodynamic  properties  of  air. 
These  equations  are 


h,  ■  /<!,, 

pl> 

(79) 

*1  *-f  <v 

pl> 

(80) 

p2  =  f  ^2  ’ 

p2) 

(81) 

u2  =  V 2  sin 

(e  -  6) 

(82) 

M2  =  V2  a2 

(83) 

al  =  /  (Ti ’ 

pi> 

(84) 

a2  =  f  (^2* 

p2) 

(85) 

Y2  =  /  ^2  ’ 

p2) 

(86) 

h  log 
Y1  ~  1  d  log 

(87) 

Inclusion  of  Equations  79  through  87  matches  the  number  of  unknowns 
with  the  number  of  applicable  equations.  A  solution  is  therefore  possible. 
The  calculation  procedure  to  be  used  is  as  follows: 


Step  1:  Find  h-| ,  p-j ,  a^ ,  and  y-j  >  using  Equations  79,  80,  84,  and  87. 

Step  2:  Assume  a  value  for  the  shock  wave  angle,  9  (this  value  must 
be  reasonably  correct  in  order  to  insure  that  the  weak  shock  solution 
will  be  determined).  The  shock  wave  angle  found  though  an  ideal-gas 
solution  will  generally  be  a  good  value  for  the  first  assumption). 


Step  3:  Find  Uj/u-j ,  using  Equation  78. 

Step  4:  Find  P^,  using  Equation  77. 

Step  5:  Find  h2»  using  Equation  75. 

Step  6:  Find  p2,  using  Equation  81. 


i 

i 
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Step  7:  Check  the  value  of  Uj/u-j  calculated  in  Step  3  with  a  value 
calculated  by  use  of  Equation  69. 

Step  8:  If  the  value  of  u^/u^  calculated  in  Step  7  is  not  close 
enough  to  the  value  calculated  in  Step  3,  assume  a  new  value  for  the 
shock  wave  angle,  and  repeat  Steps  3-7  until  the  two  values  of  Uj/u^  are 
c^se  enough.  Then  proceed  to  Step  9. 


Step 

9: 

Find 

r ^  and  aj. 

using  Equations  85  and  86. 

Step 

10: 

Fi  nd 

V 

usi  ng 

Equati on 

71  . 

Step 

11  : 

Fi  nd 

V 

using 

Equation 

70. 

Step 

12: 

Fund 

U2  ’ 

using 

the  relation:  u^  =  u^  { u^/u ^ ) 

Step 

13: 

Fi  nd 

V2’ 

using 

Equation 

82. 

Step 

14: 

Fi  nd 

m2. 

usi  ng 

Equati on 

83. 

A  subroutine  using  this  calculation  procedure  to  determine  the 
conditions  downstream  of  an  oblique  shock  wave,  given  the  upstream 
conditions  and  the  deflection  angle,  has  been  written.  The  calling 
statement  for  this  subroutine  is: 

CALL  RGSK  (DELTA,  VMA,  TA,  PA,  THETA,  VMB ,  TB,  PB) 

The  input  items  for  subroutine  RGSK  are  DELTA,  VMA,  TA,  and  PA. 

DELTA  is  the  flow  deflection  in  degrees,  VMA  is  the  upstream  Mach  number, 
TA  is  the  upstream  static  temperature  in  °R,  and  PA  is  the  upstream 
static  pressure  in  atmospheres.  THETA  is  the  shock  wave  angle  in  degrees. 
VMB,  TB,  and  PB  are  the  calculated  parameters  downstream  of  the  shock 
wave,  and  have  the  same  units  as  the  corresponding  input  parameters.  The 
information  determined  in  subroutine  RGSK  is  based  on  the  weak  shock 
solution  which  is  generally  the  solution  of  interest.  No  provision  is 
made  to  provide  the  strong  shock  solution  in  this  subroutine.  A  listing 
of  subroutine  RGSK  is  included  as  Table  A- 11.  The  use  of  subroutines 
THAIR,  NEWT,  and  CGSK  are  required  with  subroutine  RGSK. 
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3.  ISENTROPIC  EXPANSION 


The  real  gas  solution  for  a  two  dimensional  isentropic  expansion  is 
obtained  in  a  manner  similar  to  that  employed  for  the  ideal  gas  solution. 
However,  since  the  value  of  the  isentropic  exponent  (yj)  varies  throughout 
the  expansion,  the  calculation  procedure  must  be  carried  out  in  a  number 
of  smal 1  steps . 

The  applicable  equations  are 
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Using  these  equations,  and  knowing  M-j  ,  ,  and  s,  the  calculation 

procedure  is  as  follows:  • 

Step  1:  Divide  5  into  a  sufficient  number  of  equal  increments  "n". 

Step  2:  Find  Yj  =  f  (T-j  ,  P1  ) 

Step  3:  Find  v, ,  using  Equation  88 

Step  4:  Find  v.+^ ,  using  Equation  89 

Step  5:  Find  M.+,,  using  an  iterative  solution  with  Equation  88 

Step  6:  Find  T„  and  P9  ,  using  Equations  90  and  91 

l+l  i+1 

Step  7:  If  i<n  find  yt  =  (T,  )  and  repeat  Steps  2-6  until 

i  +1  ^i +1 

i=n.  The  desired  values  of  M 2 ,  ,  and  P2  have  then  been  determined. 

A  subroutine  using  this  calculation  procedure  was  written  to  calculate, 
on  a  real  gas  basis,  the  conditions  downstream  of  an  isentropic  expansion. 
The  calling  statement  for  this  subroutine  is: 

CALL  RGEXP(T1 , PI , VM1 ,DEFL,DEFLN,GAM,VMU1 ,VMU2 ,T2 ,P2  ,VM2 ) 

The  parameters  input  to  subroutine  RGEXP  are  T1 ,  PI,  VM1 ,  DEFL, 

DEF..N,  and  GAM.  If  the  calculations  are  to  be  accomplished  on  a  real  gas 
basis,  GAM  is  set  equal  to  zero.  If  ideal  gas  calculations  are  desired, 

GAM  is  set  equal  to  the  value  of  the  ratio  of  specific  heats.  T1  is  the 
upstream  temperature  in  °R,  PI  is  the  upstream  pressure  in  atmospheres, 

VM1  is  the  upstream  Mach  number,  and  DEFL  is  the  deflection  angle  in 
degrees.  DEFL  is  always  a  negative  number.  DEFLN  is  the  number  of  equal 
increments  that  the  deflection  is  to  be  divided  into  in  order  to  provide 
results  of  sufficient  accuracy.  Figure  4  indicates  the  effect  on  accuracy 
of  various  values  for  DEFLN.  This  data  was  obtained  for  a  real  gas 
case  with  VM1=6.0,  T1=2500°R,  and  PI = . 01  atmospheres.  The  total  deflection, 
DEFL,  is  15°,  and  the  resulting  values  of  M2,  T2>  and  P2  are  shown  as 
calculated  using  various  values  of  DEFLN.  For  the  Mach  6  case  shown, 
values  of  DEFLN  that  produce  incremental  deflections  less  than  about. 5° 
provide  results  that  are  generally  adequate.  For  Mach  numbers  greater  than 
6  where  the  value  of  the  isentropic  exponent  changes  more  rapidly,  a 
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2 

(R°)(10'4ATM)  M2 


DEFLN  -  Number  of  Equal  Increments 
For  Entire  Deflection 


Figure  4.  Effect  of  Calculation  Interval  on  Results  Obtained  from 
Subroutine  R6EXP 
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smaller  increment  would  be  required  for  the  same  degree  of  precision. 
Similarly,  for  lower  Mach  number  cases  a  larger  increment  would  be 
appropri ate . 

The  output  quantities  determined  by  use  of  subroutine  RGEXP  are  the 
downstream  conditions  T2,  P2 ,  and  VM2,  in  the  same  units  as  the 
corresponding  input  parameters.  Also  calculated  are  the  upstream  and 
downstream  Mach  angles  (VMU1  and  VMU2)  in  degrees.  A  listing  of 
subroutine  RGEXP  is  included  as  Table  A- 1 2 .  Subroutines  EXN  and  THAIR 
are  required  for  use  in  conjunction  with  subroutine  RGEXP. 

4.  OBLIQUE  SHOCK  WAVE  -  AX ISYMETRIC 

The  basic  Taylor-Maccoll  Theory  for  flow  over  a  cone  at  zero 
degrees  angle  of  attack  was  discussed  in  Section  II-4.  All  the  develop¬ 
ment  and  equations  through  Equation  44  are  equally  valid  for  real  gas  and 
ideal  gas  calculations.  Equation  45,  however,  is  applicable  only  to  an 
ideal  gas.  A  different  approach  will  therefore  be  required  for  real  gas 
calculations.  The  starting  point  (Equation  44)  will  be  repeated  here 
for  convenience: 


(cot  0f 


+  2u 


0 


(44) 


The  solution  is  a  continuous  function  in  the  region  between  the 
surface  and  the  shock,  so  expansion  in  a  convergent  Taylor  series  is 
possible : 


(92) 
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Using  Equations  42,  44,  92,  and  93  the  numerical  integration 
procedure,  proceding  from  the  cone  surface  (where  all  conditions  are 
specified),  to  the  shock  wave,  is  as  follows: 


Step  1 : 

Step  2: 
Step  3: 
Step  4: 


Set  =  6,  u  =  uCIIO(r,  v  -  0  (no  flow  into  or  out  of  surface), 
O  oUKr  0 

To  =  TSURF’  and  Po  =  PSURF 
Find  s,  hSURp,  pQ  =  f  (TSURFi  PSURp) 
u2 

J  ,  ,  u  SURF 

Find  ht  hSURp  +  2gJ 

Set  i  =1 


Step  5 : 

Step  6: 

Step  7: 
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Step  8: 

Step  9: 

Step  10: 

Step  1 1 : 
Step  12: 

Step  13: 

Step  14: 

Step  15: 

Step  16: 

Step  17: 

Step  18: 


Find  u.  =  u.  ,  +  (v,  ,)  + 


UeJ1 


i  -  “i_i  -  r  i  -  2 


(From  Equation  93) 


Find 


2  A  2 
u .  +  v  . 

i  l 


Fi  nd  h . 

i 


=  h 


w. 

i 


t  2gJ 


Find  T.,  Pf  =  f  (h. ,  S) 


Find  M2  =  w^  V  gyRT\ 


-1  Vi 

Find  0.  =  tan  —  +0 

i  u.  r. 

,  1  ,  1 


Fi  nd 


dp  \  _  pi-1  ~  pi 

dQ  A6 

r  j  .  r 


Find  M2  for  2-D  wedge  with  <5  =  0.  and  0=8^ 


wi 


If  M2  >  M2  ,  set  i  =  i+1  and  return  to  Step  5 


wi  1 


M?  -M? 
w.  ^i-i 

If  M,  <  M,  ,  find  0  =  0„  +  A6  .-jr - -  M  - jr— 

2w.  2i  ri-l  r  M2  ~M2  -M2  -M2. 

1  (  w1_1  w.  w.  1 


9-0r.  .  f 

Find  M2  =  M2  +  - - —  /  M  -  M? 

5  Vl  ri"  ri-l  wi  wi-l 
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Step  19:  Find  r  =  j>.  ,  +  7- - r— — 

s  l-l  5  -y_ 

ri  ri-l 
9""r.  . 

Step  20:  Find  M]  =  +  9-^-  (  ^  .  A 

t-1  r.  ri  ]  \  ^  i-l  ) 

m  computer  program  has  been  written  to  perform  the  calculation 
procedure  outlined  above.  A  listing  of  this  program  (Program  RCONE)  is 
included  as  Table  A- 1 3 .  The  use  of  subroutine  THAIR  is  required  in 
conjunction  with  Program  RCONE.  The  input  data  to  this  program  are: 

DELT,  the  cone  half  angle  in  degrees;  OTHETA,  the  calculation  interval 
in  degrees;  USRF,  the  surface  air  velocity  in  feet  per  second;  TSRF, 
the  surface  air  temperature  in  °R;  and,  PSRFA,  the  surface  air  pressure 
in  atmospheres.  The  output  for  one  set  of  input  data  is  on  a  single 
sheet.  A  sample  output  is  included  as  Figure  5.  The  input  for  this 
sample  case  is  as  follows:  DELT  =  20.,  DTHETA  =1.,  USRF  =  2000., 

TSRF  =  1000.  ,  and  PSRF  =  .5. 

A  calculation  interval  of  1°  has  been  shown  to  be  adequately  small 
for  most  ideal  gas  calculations.  (See  Figure  2.)  This  calculation 
interval  is  generally  adequate  for  real  gas  calculations  as  well.  If 
extreme  precision  is  required,  a  smaller  interval  can  be  selected. 

The  calculation  procedure  and  computer  program  discussed  in  the 
preceding  paragraphs  is  an  inverse  approach  to  the  problem.  The  condi¬ 
tions  on  the  cone  surface  have  been  specified  and  the  conditions  upstream 
of  the  shock  wave  are  calculated.  In  most  cases,  however,  the  data 
available  is  the  conditions  upstream  of  the  shock  wave,  and  the  informa¬ 
tion  to  be  determined  is  the  conditions  downstream  of  the  shock  wave. 

It  is  possible  to  obtain  the  desired  information  through  iterating  with 
program  RCONE.  In  some  cases  this  may  be  the  most  appropriate  technique. 
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However,  a  computer  subroutine  has  been  written  to  perform  the  iterative 
computation  procedure  required  with  the  generally  available  data.  The 
calling  statement  for  this  subroutine  is: 

CALL  RGCON  (VM1,  T1 ,  P1A,  CELT,  THE,  THETAS,  DSR ,  VM2C ,  PHIC,  PRC,  TRC , 

VM2S ,  PRSHK,  TRSHK,  PHIS,  VM2SRF,  PRSRF,  TRSRF ) 

The  first  five  items  in  this  sequence  are  inputs.  VM1  is  the 
upstream  Mach  number,  T1  is  the  upstream  static  temperature  in  °R,  P1A 
is  the  upstream  static  pressure  in  atmospheres,  and  DELT  is  the  cone  half 
angle  in  degrees.  The  fifth  input,  THE,  has  been  included  to  provide 
data  on  one  specific  ray  lying  between  the  cone  surface  and  the  shock 
wave.  THE  is  the  ray  angle,  in  degrees,  for  which  data  is  required.  The 
parameters  found  by  use  of  subroutine  RGCON  are  as  follows:  THETAS  is 
the  shock  wave  angle,  DSR  is  the  non-dimensi onal i zed  entropy  increase 
across  the  shock  wave,  VM2C  is  the  Mach  number  on  the  specified  ray, 

PHIC  is  the  flow  angle  on  the  specified  ray,  PRC  is  the  ratio  of  the 
static  pressure  on  the  specified  ray  to  the  upstream  static  pressure, 

TRC  is  the  ratio  of  the  static  temperature  on  the  specified  ray  to  the 
upstream  static  temperature,  VM2S  is  the  Mach  number  just  downstream 
of  the  shock  wave,  PRSHK  is  the  static  pressure  ratio  across  the  shock 
wave,  TRSHK  is  the  static  temperature  ratio  across  the  shock  wave,  PHIS 
is  the  flow  angle  just  downstream  of  the  shock  wave,  VM2SRF  is  the  Mach 
number  on  the  cone  surface,  PRSRF  is  the  ratio  of  the  static  pressure 
on  the  surface  to  the  upstream  static  pressure,  and  TRSRF  is  the  ratio 
of  the  static  temperature  on  the  surface  to  the  upstream  static  temperature. 
All  angles  are  in  degrees.  A  listing  of  subroutine  RGCON  is  included 
as  Table  A- 1 4 .  The  use  of  subroutines  THAIR  and  NEWT  is  required  in  con¬ 
junction  with  subroutine  RGCON. 

To  illustrate  the  use  of  subroutine  RGCON,  and  additional  program 
has  been  written.  The  listing  of  this  program,  which  is  entitled 
program  RGCONE,  is  included  as  Table  A- 1 5 .  The  input  data  to  this 
program  are:  VM1 ,  the  upstream  Mach  number;  T1 ,  the  upstream  static 
temperature  in  °R;  P1A,  the  upstream  static  pressure  in  atmospheres; 

DELT,  the  cone  half-angle  in  degrees;  and  THE,  the  angle  in  degrees  of 
a  ray  for  which  data  is  desired.  A  sample  output  for  program  RGCONE  is 


AFAPL-TR-79-2023 


included  as  Figure  6.  The  input  data  used  to  generate  Figure  6  is  as 
follows:  VM1 =3.0 .  T1=390°R,  P1A=1 .0  atm,  DELT=20 ,  THE=22.  This  data 
is  similar  to  that  used  by  program  CONES  with  subroutine  CONSK  to  obtain 
the  printout  presented  in  Figure  3  for  ideal  gas  with  ratio  of  specific 
heats  equal  1.4.  It  can  be  seen  by  comparison  of  Figures  3  and  6 
that  the  results  are  very  much  in  agreement.  This  is  to  be  expected  for 
the  example  shown.  Note,  however,  that  the  computer  time  required  to 
generate  Figure  6  was  in  excess  of  50  seconds,  while  the  computer  time 
for  Figure  3-type  calculations  is  nearly  insignificant. 
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Figure  6.  Sample  Output  for  Program  RGCONE 


AFAPL-TR-79-2023 


SECTION  IV 

EXAMPLES  USING  DEVELOPED  TECHNIQUES  AND  SUBROUTINES 

The  utility  of  the  various  computation  techniques  and  computer 
subroutines  discussed  in  the  preceding  sections  will  be  demonstrated 
in  this  section.  Three  typical  supersonic  airflow  problems  will  be 
employed  to  accomplish  the  demonstration.  Each  problem  will  be  approached 
in  three  ways.  First,  for  a  low  Mach  number  case,  hand  calculations, 
using  the  charts  in  Appendix  B,  will  be  described.  Second,  a  computer 
program,  written  to  use  the  ideal  gas  subroutines  of  Section  II,  will  be 
described  and  demonstrated  for  low,  intermediate,  and  high  Mach  number 
cases.  Finally,  a  computer  program,  written  to  use  the  real  gas  sub¬ 
routines  of  Section  III,  will  be  described  and  demonstrated  for  the  same 
Mach  numbers.  A  comparison  of  the  final  two  techniques  will  illustrate 
the  validity  of  the  ideal  gas  techniques. 

Note  that  all  the  techniques  discussed  in  this  document  are  for 
inviscid  air  flows  only.  Inclusion  of  viscous  effects  and  correction 
techniques  is  not  within  the  planned  scope  of  this  document.  In  reality, 
all  airflows  are  dependent,  to  a  greater  or  lesser  degree,  on  viscous 
effects.  For  this  reason,  any  comparisons  of  results  using  the  tech¬ 
niques  discussed  in  this  document  with  actual  test  data  are  not  valid. 

The  first  of  the  supersonic  airflow  problems  to  be  considered  is  a 
two-dimensional  swept  wing  of  infinite  span.  The  procedure  for  deter¬ 
mining  the  lift  and  drag  of  such  a  wing  will  be  discussed.  The  second 
problem  concerns  a  two  dimensional  supersonic  inlet  for  an  airbreathing 
propulsion  device.  The  back  pressure  required  to  maintain  the  normal 
shock  on  the  inlet  lip, and  the  total  pressure  recovery  just  downstream 
of  the  normal  shock  will  be  calculated  for  an  inlet  with  a  two-step  ramp, 
at  Mach  numbers  on  or  below  the  design  Mach  number.  The  third  problem  is 
similar  to  the  second,  but  uses  a  conical  centerbody  inlet  at  zero  angle 
of  attack.  Again,  the  back  pressure  required  to  keep  the  normal  shock  on 
the  inlet  lip  and  the  total  pressure  recovery  just  downstream  of  the  normal 
shock  will  be  calculated  for  cases  at  or  below  the  design  Mach  number. 
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1.  PROBLEM  1  -  DETERMINATION  OF  PRESSURE  DISTRIBUTION  AND  LIFT 
AND  DRAG  FOR  A  SWEPT  WING  OF  INFINITE  SPAN 

The  basic  features  and  dimensions  of  the  wing  being  considered  are 
illustrated  in  the  following  diagrams: 


wing  cross  section 
parallel  to  flight  direction 


wing  viewed 
from  above 


The  parameters  to  be  specified  for  this  problem  are:  the  upstream 
Mach  number,  static  temperature,  and  static  pressure,  MQ,  Tq,  and  Pq; 
the  wing  angle  of  attack,  or,  and  the  wing  dimensional  quantities,  y^,  y2» 
°1 ’  a2’  a3’  a4’  and  The  Pressure  distribution  on  the  wing  upper 
and  lower  surfaces,  and  the  wing  lift  and  drag  are  to  be  calculated. 

For  the  case  demonstrating  the  hand  calculation  technique,  the 
specifications  will  be 

M  =  2.5,  T  =  500°R ,  P  =  .5  atm,  a  =  2° 

0  0  0 

y1  =  .2  ft,  y2  =  .5  ft 

o]  =  5°,  a2  =  15°,  a3  =  20°,  a4  =  10° 

A  =  60° 

The  first  step  in  the  calculation  procedure  is  to  determine  the 
equivalent  upstream  conditions  due  to  the  wing  sweep.  From  Figure  B-9, 
we  find 


en 

jjj-2-  =  .867  -  Me  =  2.5  (.867)  =  2.168 
0  0 
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And,  from  Figure  B- 1 0 


a  =  2.3° 

e 


a  =  17.2° 
e2 


a  =  22.8° 

e3 

a  =  11.5° 

e4 

To  determine  the  pressures  and  forces  on  the  wing  surfaces,  one  must 
first  determine  the  deflections  the  airflow  will  encounter.  The  deflections 
and  forces  are  shown  in  the  following  diagram: 


wing  cross  section 
normal  to  leading  edge 
The  deflections  then  are 


5  =  a  -  a  =  5.7°  -  2.3°  =  3.4° 

el  el  e 

<5_  =  a  +  a  «  17.2°  +  2.3°  =  19.5° 

e2  e2  e 

=  -  cr  -  a  =  -  5.7°  -  22.8°  =  -28.5° 
e3  el  e3 

6  =  -  -  o  =  -  17.2°  -  11.5°  =  -  28.7° 

e4  e2  e4 

The  next  step  is  to  determine  the  unspecified  dimensions  of  the  wing. 

Referring  to  the  diagram  at  the  beginning  of  Section  IV. 1: 


x 


1 


^1 

tan 


■  b 

tan  5° 


=  2.286  ft 
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y2  5 

=  ■  -  -  -  ■  =  v-  =  1.866  ft 

tan  tan  15 

y3 

tan  a3 

(94) 

y4 

tan  o4 

(95) 

+  y4  =  yl  +  y2 

(96) 

+  x3  =  x2  +  x4 

(97) 

Equations  94  through  97  can  be  solved  for  to  find: 

yl  +  ^2  +  ^x2  ‘  X1 ^  tan  °4 
x3  tan  a 3  +  tan  c?4 

=  .2  +  .5  +  (1  .866  -  2.286)  tan  10°  =  158 

3  tan  20°  +  tan  10° 


Then,  using  Equations  97,  95,  and  96: 

x4  =  2.286  +  1 .158  -  1 .866  =  1 .578  ft 

y4  =  1 .578  tan  10°  =  .278  ft 

y3  =  .2  +  .5  -  .278  =  .422  ft 

The  total  cord  of  the  wing  is: 

X  =  x,  +  x3  =  2.286  +  1.158  =  3.444  ft 

To  calculate  the  pressure  distribution  over  the  wing,  a  sketch 
indicating  the  approximate  location  of  shock  waves  and  expansion  fans 
will  be  hel pful . 
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Since  we  are  interested  only  in  the  forces  on  the  wing  surfaces 

there  is  no  need  to  calculate  the  oblique  shock  wave  angles  or  the  Mach 

line  angles  for  the  expansions.  In  addition,  there  is  no  reason  to 

consider  the  trailing  edge  oblique  shock  waves.  The  Mach  number  change 

and  static  pressure  ratio  across  each  of  the  leading  edge  shock  waves 

can  be  determined  by  use  of  Figures  B-4  and  B-6.  However,  we  should 

first  check  to  make  sure  that  the  deflections  are  not  so  great  as  to 

preclude  an  attached  shock  wave.  This  is  done  with  Figure  B-2.  Using 

Figure  B-2,  we  find  that  at  M  =  2.168  the  maximum  deflection  is  25.6°. 

e 

o 

Both  of  the  leading  edge  deflections  are  less  than  25.6°,  so  the  oblique 
shock  waves  will  be  attached.  Then,  using  Figures  B-4  and  B-6,  we  find: 


And: 


.645 


M1  =  .94  (2.168)  =  2.038 
M2  =  .645  (2.168)  =  1 .398 

P]  =  ,81 9  =  a^m 
P2  =  if 9  =  1 -^33  atm 
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The  expansions  at  the  top  surface  and  bottom  surface  corners  can  be 
addressed  by  use  of  Figure  B- 11 .  The  Prandtl -Meyer  angles  upstream  of  the 
corners  are 


'V1 

v2 


27.3° 

8.9° 


Downstream  of  the  corners. 


v 


3 


v4  =  v2  " 


27.3°  +  28.5°  =  55.8° 
8.9°  +  28.7°  =  37.6° 


Again  using  Figure  B-ll, 

M  =  3.33 
3 

M4  =  2.44 

And,  since  the  expansion  flow  is  isentropic. 


P 


3 


+ 


+ 


For  y  =  1.4: 


1 


1 


+ 


+ 


P 


3 


1  +  .2  (2.038)2 
1  +  .2  (3.33)2 


3.5 


.085  atm. 


P 


4 


1  .433 


1  +  .2  (1.398)2 
1  +  .2  (2. 44) 2 


3.5 


.292  atm. 
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All  the  forces  on  the  wing  surfaces  have  now  been  defined.  Ail  that 
remains  is  to  sum  these  forces  in  the  appropriate  orientation.  The  force 
normal  to  the  wing  as  illustrated  by  the  figure  on  page  52  is  defined: 

N  =  P2  x2  +  p4  x4  -  P1  x1  -  P3  x3  (98) 

The  force  along  the  wing  axis  is  defined: 

F  =  P3  y3  +  P4  y4  ‘  P1  yl  '  P2  y2  (99) 

Using  Equations  98  and  99, 

/ 

N  =  21 1 6 [ 1 .433  (1.866)  +  .292  (1.578)  -  .641  (2.286)  -  .085  (1.158)] 

=  3469.3  lbs/ft  of  wing  span 

F  =  2116  [.085  (.422)  +  .292  (.278)  -  .611  (.2)  -  1.433  (.5)] 

=  -1527.0  lbs/ft  of  wing  span 

The  lift  and  drag  can  be  found  from  the  normal  and  wing  axis  drag  by 
using  the  equations: 


Lift  =  N  cos  a  +  F  sin  a  (100) 

Drag  =  N  sin  a  -  F  cos  a  (101) 

Then:  Lift  =  3469.3  cos  2°  -  1527.0  sin  2°  =  3413.9  lbs/ft  of  wing  span 

Drag  =  3469.3  sin  2°  +  1527.0  cos  2°  =  1647.1  lbs/ft  of  wing  span 

A  computer  program  designed  to  accomplish  the  calculations  described 
above  has  been  written.  A  listing  of  this  program  is  included  as  Table  A- 1 6 . 
The  input  for  this  program  is:  MO,  TO,  and  PO,  the  Mach  number,  static 
temperature,  and  static  pressure  upstream  of  the  wing:  Yl  and  Y2,  the 
thickness  of  the  wing  above  and  below  the  centerline;  SIG1 ,  SIG2,  SIG3, 
and  SIG4,  the  angle  of  each  wing  surface  with  respect  to  the  wing  centerline; 
LAMBDA,  the  wing  sweep  angle;  and  ALPHA,  the  wing  angle  of  attack. 

Temperature  is  in  °R,  pressure  is  in  atmospheres,  dimensions  are  in  feet, 
and  angles  are  in  degrees.  The  output  of  this  program  is  a  single  sheet 
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for  each  input  case.  Output  for  three  sampl  e probl ems  is  included  as 
Figures  7,  8,  and  9.  The  first  case  (Figure  7)  has  input  data  identical  to 
the  preceding  hand  calculated  example.  The  second  case  as  illustrated 
in  Figure  8  is  for  a  wing  at  a  Mach  number  of  6.0.  The  final  case  is  for 
a  wing  at  a  Mach  number  of  12.0.  This  case  is  illustrated  in  Figure  9. 

The  Mach  6.0  and  Mach  12.0  cases  illustrated  are  actually  not  valid 
since  the  calculations  have  been  made  on  an  ideal  gas,  y  =  1.4,  basis. 

The  degree  of  invalidity  will  be  demonstrated  by  comparison  with  results 
obtained  on  a  real  gas  basis  for  identical  input  data.  In  order  to 
accomplish  this  comparison,  the  swept  wing  computer  program  has  been 
modified  to  include  real  gas  subroutines  RGSK  and  RGEXP.  A  listing  of 
the  revised  computer  program  is  included  as  Table  A- 1 7 .  The  same  input 
sets  used  to  generate  Figures  7,  8,  and  9  have  been  used  with  the  revised 
program  to  generate  Figures  10,  11,  and  12.  Comparison  of  the  correspond¬ 
ing  output  sheets  indicates  that  for  the  Mach  2.5  example  the  results 
are  very  similar.  At  Mach  6.0  there  is  some  difference,  and  at  Mach  12.0 
there  are  major  discrepancies. 

2.  PROBLEM  2  -  DETERMINATION  OF  BACK  PRESSURE  AND  PRESSURE 

RECOVERY  FOR  A  TWO-DIMENSIONAL  INLET 

The  inlet  to  be  considered  is  as  described  in  the  following  sketch: 


T 

o. 


al 

To  avoid  undue  complication  a  number  of  restrictions  will  be  placed 
on  the  cases  to  be  addressed.  The  oblique  shock  w «ves  will  not  intersect 
within  the  capture  area  of  the  inlet,  and  neither  of  the  oblique  shock 
waves  will  enter  the  inlet;  the  angle  of  attack  will  be  less  than  the 
first  ramp  angle;  neither  of  the  ramp  angles  will  be  so  large  that 
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Figure  7.  Mach  2.5  Case  for  Ideal  Gas  Swept  Wing  Problem 
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Figure  8.  Mach  6  Case  for  Ideal  Gas  Swept  Wing  Problem 
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Figure  9.  Mach  12  Case  for  Ideal  Gas  Swept  Wing  Problem 
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Figure  10.  Mach  2.5  Case  for  Real  Gas  Swept  Wing  Problem 
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Figure  11.  Mach  6  Case  for  Real  Gas  Swept  Wing  Problem 
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attached  shock  waves  are  precluded;  and,  the  inlet  back  pressure  will  be 
such  that  the  normal  shock  is  maintained  exactly  at  the  cowl  lip.  The 
assumption  of  inviscid  flow  will,  of  course,  also  apply  to  this  problem. 

The  parameters  to  be  specified  for  this  problem  are:  the  upstream 

Mach  number,  static  temperature,  and  static  pressure,  M  ,  T  ,  and  PQ; 

the  inlet  angle  of  attack,  a;  and  the  ramp  angles,  and  a The 

static  pressure  at  station  3  (inlet  back  pressure),  and  the  inlet  total 

pressure  recovery  (P.  /P  )  are  to  be  calculated. 

r3  ro 

For  the  case  demonstrating  the  hand  calculation  technique,  the 
specifications  will  be 

Mq  =  2.5,  Tq  =  500°R ,  Pq  =  .5  atm 

o i  =  10°,  a2  =  20° 

a.  =  -5° 

The  flow  deflections  over  the  ramp  are  easily  determined  as  follows 

*1  =  °1  -  *  =  1°°  +  5°  "  150 
62  =  °2  "  °1  =  20°  '  10°  =  10° 


Using  Figures  B-3,  B-4,  B-6,  and  B-7,  we  find 


<5-,  M,  Asn  1 

=  -406,  —  =  .749,  ^  =  .405,  — =  .0738 

91  =  =  36-95°  _5°  =  31 -95° 

M]  =  .749(2.5)  =  1 .873 

^1  =  405  =  1  ,333  atm 


Again  using  the  Appendix  B  figures,  we  find: 


62  M 

—  =  ?37  — 
92  M] 


.812, 


.597, 


As 


1-2 


.0145 
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And 


10° 

.237 

.812 

1  .235 
.597 


+  oj  =  52.19° 
(1.873)  =  1 .521 
=  2.069  atm. 


Finally,  Figure  8-1  is  used  to  determine  the  conditions  following 
the  normal  shock.  From  Figure  B-l , 


And 


.395, 


As2-3 

R 


.0798 


„  2.069  r 

^3  =  — 395  5.238  atm. 


This  is  the  back  pressure  (crea^d  by  combustion  of  fuel  in  the 
engine)  that  is  required  to  maintain  the  normal  shock  wave  exactly  at 
the  inlet  lip.  A  higher  back  pressure  would  force  the  normal  shock 
forward  of  the  lip  and  cause  spillage  of  air.  A  lower  back  pressure  would 
allow  the  normal  shock  to  move  into  the  inlet  and  a  complicated  pattern 
of  reflected  shocks  and  expansion  waves  would  then  precede  the  normal 
shock. 


The  dimensionless  entropy  increase  for  the  inlet  is  determined  by 
summing  the  increases  for  each  of  the  shock  waves. 


As0-3  As 0-1  ,  As1-2  .  As2-3 

R  =  R  +  R  +  R 


.0738  +  .0145  +  .0798  «  .1681 


The  pressure  recovery  of  the  inlet  is  easily  found  from  the  rela¬ 
tionship 


-.1681 

e 


.8453 
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A  computer  program  has  been  written  to  accomplish  the  calculations 
described  above.  A  listing  of  this  program  is  included  as  Table  A- 1 8 . 

The  inputs  for  this  program  are  MO,  TO,  PO,  the  Mach  number,  static 
temperature  in  °R,  and  static  pressure  in  atmospheres,  upstream  of  the 
inlet;  SIG1  and  SIG2,  the  two  inlet  ramp  angles,  in  degrees;  and,  ALPHA, 
the  inlet  angle  of  attack  in  degrees.  The  output  of  this  program  is  on 
a  single  sheet,  and  includes  both  the  ideal  gas  and  the  real  gas  solutions 
appropriate  to  the  input  data.  Output  for  three  sample  cases  is  included 
as  Figures  13,  14,  and  15.  The  first  case  (Figure  13)  has  input  data 
identical  to  the  preceding  hand  calculated  example.  The  second  case,  as 
illustrated  in  Figure  14,  is  for  an  inlet  at  a  Mach  number  of  6.0.  The 
third  case  is  for  an  inlet  at  Mach  12.0.  The  output  sheet  for  this  case 
is  Figure  15.  The  degree  of  discrepancy  between  real  gas  and  ideal  gas 
calculations  can  be  observed  at  each  of  the  three  Mach  numbers. 

3.  PROBLEM  3  -  DETERMINATION  OF  BACK  PRESSURE  AND  PRESSURE 

RECOVERY  FOR  A  CONICAL  INLET 

The  inlet  to  be  considered  for  this  problem  is  described  in  the 
following  sketch: 


The  restrictions  on  this  problem  include  the  following: 

(1)  The  inlet  is  at  a  zero  angle  of  attack. 

(2)  The  shock  wave  angle  (0)  is  greater  than  the  cowl  angle 

(0COWL^ 

(3)  The  cone  half  angle  (6)  is  small  enough  that  an  attached 
shock  wave  is  possible. 

(4)  The  back  pressure  is  regulated  so  that  the  normal  shock  wave 
is  exactly  on  the  cowl  lip. 
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BACK  pressure  AND  PRESSURE  RECOVERY 
POP  2-D  TNLCT 
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Figure  13.  Mach  2.5  Case  for  2-D  Inlet  Problem 
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Figure  14.  Mach  6  Case  for  2-D  Inlet  Problem 
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Figure  15.  Mach  12  Case  for  2-D  Inlet  Problem 
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(5)  The  properties  of  the  air  upstream  of  the  normal  shock  can  be 
adequately  represented  by  the  conditions  on  a  ray  midway  between  the  cone 
surface  and  the  cowl . 

(6)  The  airflow  is  inviscid. 


The  specified  data  for  this  problem  are:  the  upstream  Mach  number, 
c  temperature,  and  static  pr 
angle,  5;  and,  the  cowl  angle,  0, 


static  temperature,  and  static  pressure,  M  ,  Tq  ,  and  PQ;  the  cone  half 

The  static  pressure  at  station  2 


COWL' 

(inlet  back  pressure)  and  the  inlet  total  pressure  recovery  (P.  /P.  ) 
are  to  be  calculated. 


i 


For  the  case  demonstrating  the  hand  calculation  technique,  the 
specifications  will  be 


Mo  =  2.5,  To  =  500°R ,  PQ 
6  =  15°,  eC0WL  =  25° 


.5  atm 


The  shock  wave  angle  is  found  by  use  of  Figure  B-13  as  follows: 
|  =  .527  -  0  =  -~y  =  28.46° 


Using  Figure  B-14,  we  find 


And 


0r=15° 


.847  -+  M,  =  2.118 

V15° 


9  =28.46° 
r 


.896  -*■  M,  =  2.240 

6  =28.46° 
r 


Since  we  cannot  find  the  properties  on  specific  rays  between  the  k.one 
surface  and  the  shock  wave  by  use  of  the  available  figures,  the  further 
assumption  will  be  made  for  the  hand  calculations  that  the  properties 
vary  linearly  from  the  cone  surface  to  the  shock  wave.  The  Mach  number 
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on  the  ray  midway  between  the  surface  and  the  cowl  (-3  =20°)  is  then 

rMP 

found  as  fol lows : 


M 


1„ 

u 


M1  '  M1 

j  =28.46”  u  =15° 

r  r 


Ml  =  2.118  +  28^46~-TF  ^2 -24°-2 • 1 1 8)  =  2-163 

8  >- 
rMP 


To  find  the  static  pressure  downstream  of  the  shock  wave,  we  use 
Figure  B-16: 

Ip. 


1 


=  .555  P 


.5 


Qr=15° 


1  555 

'e  =15°  ‘ 

r 


=  .9009  atm. 


=  .673  -*  P 


.5 


6r=28.46' 


1 8  =28.46°  ,673 

r 


=  .7429  atm. 


Again  using  the  linear  variation  assumption: 


P,  =  .9009  +  (.7429  -  .9009)  =  .8422  atm. 

V  ,3-16 

rMP 


The  dimensionless  entropy  increase  is  found  using  Figure  B-17: 


.007 


Using  properties  on  the  ray  midway  between  the  surface  and  the  cowl 
as  upstream  conditions,  the  conditions  downstream  of  the  normal  shock  are 
found  using  Figure  B-l : 

M2  =  .552 


P 


1 


.1894  -  P2 


.8422 

.1894 


=  4.447  atm. 
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.4370 


The  value  of  P2  determined  above  is  the  back  pressure  required  to 
maintain  the  normal  shock  wave  exactly  at  the  inlet  lip. 

The  overall  entropy  increase  is  found  by  summing  the  losses  through 
the  conical  shock  wave  and  the  normal  shock: 


AsQ-2  =  As0-1  As1-2 

R  "  R  R 


.007  +  .437  =  .444 


The  pressure  recovery  then  is: 


.6415 


A  computer  program  has  been  written  to  accomplish  the  calculations 
described  above.  The  only  deviation  from  the  illustrated  procedure  occu 
because  it  is  possible,  using  the  cone  subroutines,  to  find  the  conditio 
on  the  specific  ray  that  falls  midway  between  the  cone  surface  and  the 
cowl.  A  listing  of  the  computer  program  is  included  as  Table  A-19.  The 
required  input  data  for  this  program  are  MO,  TO,  and  P0,  the  conditions 
upstream  of  the  inlet;  BELT,  the  cone  half  angle  in  degrees;  and  CWLANG 
the  cowl  angle  in  degrees.  The  output  of  the  program  is  a  single  sheet 
and  includes  both  the  ideal  gas  and  the  real  gas  solutions  appropriate 
to  the  input  data.  Output  for  three  sample  cases  is  included  as  Figures 
16,  17,  and  18.  The  input  data  for  the  first  case  (Figure  16)  is  identi 
to  the  conditions  for  the  hand  calculated  example.  The  second  case,  as 
illustrated  in  Figure  17,  is  for  an  inlet  at  a  Mach  number  of  6.0.  The 
third  case  (Figure  18)  is  for  an  inlet  at  Mach  12.0. 
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• ruL  GAS  SOLUTION 

11=  2.14  Pl  =  *>.  73  7^ -C1AT<  71=  5f  7,5h 

12=  .5o  P2  =  4.534f*"GA-t*  ’  ?=  10  51.  5R 

THE"  AS  =  2«.  4P3EG 

'  lQ:J I '  E'J  2Au<  nkESbUh  L  =4.534l  ♦  00  A  TK 
-’RESSJRE  kEGOVEFY  (FT  27  PT  0)  =  .6452 

Figure  16.  Mach  2.5  Case  7ur  Conical  Inlet  Problem 
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U=  4.3’  31  =  2.  3  35-.-niG  .  H  11*  bbl.2r 

nz-  .43  3  2*5. 2  32-*v  C A 1 1  T2  =  3uSb.0R 

T HETLG=1  ). 0  LCEG 

\-0'Jl  PpCK  J<- ->SUEE=5.  232E+C1  C  Al  M 
36  E  -SURE  f  EGO  V£<Y(PT  2/P  ft)  =.l<  7  51 

i  t  al  g  .  >  :oLJtio*. 

il=  4.  3 3  °J  =2. 3  55£-0  ia  ;  H  .1=  bfcl.21 

12=  .41  3  2=5  •  35  }E+  '  (I*  i  M  1,a  2043.4P 

'h  ■  :  t  i=i n.  j  ir e  G 

<£QOI  El’  3  ACK  P f  tG SURE  =5 . 350  t*  C  C  A  .  M 

3lvE  .SURE  F  EGO  V£  ?  Y  (P  f  2/PT  t)  *.  0  63b 

Figure  17.  Mach  6  Case  for  Conical  Inlet  Problem 
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2 a : <  pi  isi l'--.  an j  f  ■* -. .1  u  i  ^c’v  rv 

.*v»  CU;-ii  CAL  r  *  >.  c.  < 

s  12  3  P3-- 5-  c;t>,3A  M  J=  9G5*  £  t 

rSLTAsiS.TDEC 

CONL  A  N  G  L  *•  *  ll.CO'.f  C> 

TO*:Al  jAv  L  3 J  -  Iv>S 

..i=  7.99  F  1-3.7  8  4  ft  2A  Th-  Tl=  1(81,  If. 

M2-  .39  F2f  2. 6  1  5E*90  ATH  7  2*  l^t  ‘  54 . 9K 

“HFTA2-1I. 90LCG - '  ~  - 

REQUIRED  3ACK  P RPSSURt-  2.  o  1 5Z  +  J  -J A  TM 

P'CCHK:  rECOVE/Y (P~2/PVu )-  0C47 

A  L  > > A  i  L3LU‘.  i  ON 

Mi=  7.99  Pl  =  3.784f;2ATM  Tl  =  1.81.1k 

”2 -  .32  P2^S.iT2tre*ulU7*  ~  r?*  rrrrj-srr 

-HLTAS=11. 8SLEG  - - 

Pf-QUlREtTEACK  P  RE  SSURCsT.X  ZZZ+V  ?ATW  ~ 

PTEVSURE  PTCOVEf  t  (P7-2/PT0T=rTnri  3- - 

Figure  18.  Mach  12  Case  for  Conical  Inlet  Problem 
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APPENDIX  A 

COMPUTER  PROGRAM  LISTINGS 

Listings  of  computer  programs  and  subroutines  referred  to  in  the  text 
are  included  in  this  appendix. 
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TABLE  A-l 

SUBROUTINE  FOR  NORMAL  SHOCK  CALCULATIONS 
(IDEAL  GAS) 


J*Z  -  S3  ( \  CoA"-i  . )  2*2  . )  /  rt*  <v,AH-i  . )  ) ) 

I  2  -  -1**  t-^i- )♦(  l  GAr,  +  1 . ) 

j.  *  *  L 

r1  L  -  :  i»  (c  *2-  jlf.*- 1 .  )  /  U»AM+1.  ) 

c.  i  '.*% 

£  N  j 
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TABLE  A-2 

SUBROUTINE  FOR  OBLIQUE  SHOCK  CALCULATIONS 
(IDEAL  GAS) 


0  J9'  .nil  .  wOjK*  l  ,  :  L  ilfi  ijf  Ijjm  1»3*(*_lT  'jlri-T  (/f1  tjIujPc) 

(jfcL  ■ 

<  =  i..-c  t) 

V  =  l/'l-MMoii-l.lfc. 

2  =  ;  "  c)  *  .  )  +2. 

jj  -  _  ?  . '  (Zt*2)/{-J1_.*Xi+s3;>.’'Y*Z/(‘»3-.MX*»2>)-/’2^.*ir**2>/l-,32.* 
1  (<**)  ) 

j  L  -  _  ?  .  *  Y *(2**3)  /  l-. it.  (  X  +  n  ;  ) 

Tv.l*  =  .’aOI/Z.) 

JILT  -  -  i  «  M  i  i  /  3  J  *  j  ■  )  /  ,  )  1  ^  -t  „  o 

1  r  (  *  )  1 1  2  i  2 

'  I  H  l  1  ~  =  '  . 

J  1  C.  -  si  . 

T  _  =  . 

D  2  -  • 

~  '.iliv  % 

i  I  r  r s.  l  r  - . ;  i )  i ,  3 ,  • 

■  T  n  ■_  T  i  -  2  >  An  1 3  i/jjm  .**3-i.J)/.')17-c>33 

\I"Z  =  Yii 
T  t  =  ’  i 
JZ  =  ->1 

-s  t”  i  u  ‘  • 

U  0 1  -i  -  C  C  j  ( L  l  L  )  /  j  1  N  (  3  f  l  I 
t'  =  ? .  * :  o*  x/y 
i  - 

*  -  2  ,*Cc.D/y 
«  =  (  -,*U-F**2)  /  3  . 

3  =  (l  .*  (I  ♦*  3) -9.*,=  *'j+2  .*0/2?. 

IF  (ri)  3»  j|t 

JoPH  -  T  < -2/ .  *(  o*' 2)  /  <*  ,  4  <***3) )  ) 

Ml  =  an  IV  (  IjCKTU.  ) /LOrH) 

G  ^  VT  7 

f  „UPH  =  -^ur.I  l-?7.*  /  (  ♦.*  (A**0>)  ) 

PH  -  •  .1  rlT-9-aTAMt  (  i  l.-.JLr  rl**2))  A -CiPH)  ) 
isF  (  1  ~  -  1  )  £  I  Q  )  9 

is  X  =  ,.»£^T(-H/j.)*iiO:(^.i3o.MPn/3.i  ) 

5  s3  “  J  1  .• 

3  EX  =  £.*_s«r\1  { - A/ 3 , ) *  COS  * ^ h/ 3  •  ) 

1  1  rlr  i  =  ATaN(tX-P/3  ,  ) 

7Hci>  =  i  hi_T/  .  f  1/  *-  y  j  3 
T  r;b  =  U-  A*SiN(THE  i  )  )♦*. 

V-  2  =  Si-.T  ( (  (C-A  1-1,  )  *T‘L*L,  )  /  i  (  (Sltf  <  Ti ET-Jtu ) ) **2> ♦  U.’GAf1  *  i  MS-GA.i 
1» 1.  >  )  ) 

12  =  .  i*  (  U  .*  GArtM  ‘1S-GA"'**.)  *  (  (GAi-i.  P  TM3*2.  )  )  /  U  (GAM+1 .)  **2>  *T  MS  ) 
°2  =  -l*  <  <2.*Hl*T  ib-U  .♦i.j /(GArt  +  i.  n 
*<tru  ?.m 
END 
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TABLE  A- 3 

SUBROUTINE  FOR  SWEPT-WING  CONVERSIONS 


SUBROUTINE.  SWLP(AwAH,VM,*L:>--l-i  »il  &  » VMEU*  ALeE^»SIGEQ) 

p  n-  niTt1 

VhU  =  ASIN(l./Vr.)/RD 
IF<V^U.oT , Ml AM) SJ  TO  1 

VMt Q=  VM*  SQRT  ( 1 .  - (  C  JS  ( A. °HA *  \ J )  *COS(AlAM*!R-j)  )**2) 

AL?  E0  =  ATAM3  IN  (  AlFHa*RO)  /  <  3 IH  <  A  Lm  M*kD )  *  GOS  UL  PH  A*  R  J )  >  )  /  RG 
aIGEQ  =  ATAMSIN(SI3*.<.3)/  iSil  (At-AM  +  kQ)  *CQ S < Si G*kO ) >  )/KD 
GO  TO  2 

1  HRI1F(6,T)VMU,ALA« 

3  FORMAT ( 1H  ,‘hACH  hNIGLL  jREAIEP  THA,<  SWEEP  AH&lE  ( <iU  =  *  » 1P£  1 J  •  2  ,  b  X ,  v 
lLAHBnA=*f E1P. £,♦)•> 

2  RETURN 
ENO 
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TABLE  A- 4 

SUBROUTINE  FOR  ISENTROPIC  TURN  CALCULATIONS 
(IDEAL  GAS) 

sueFDUTi^t  -xf-  (j- 
c<Ti9,*(«L  <  x\‘C 
con-u./?* !:C1/  &i,  j? 

•V_  =  .  Jl  7  UR  ” 

=  L>r  t. *r  -■ 

cl  -  Sj~"  US*  >"H.)'CC, 

v/'jl  -  SI  *  i'  i  AN<  S  3  <  r  (  V*  1»*  2-1.  1  /SI ) -AT  AN<  i  (  V  '1*  ‘2-  1.  )  ) 
V’-J:  =  V '!  'J  1  *  r  F  L 

V/'-  ?',=  ’  ul 

j--, x-mn  i . 

«  - 1 . 

".CS V  =  1.  E-  c 

0rLT"=»31 

3  A  L  .  iE  W'  i(V*ic,  VI  •!,  7  IX  ,  CCMJ  i  r>  3  ,CrLl'4 , 4H  TXi'l''  »uXNC  ,  X) 

T2=il‘(l.<-(G^'-l.)*(v/‘'l**2>'2.)/(l.  +  (GAl-l.)*W<2**2)/2.) 
Pi  -  i*  r  2/f  1)  »♦  (  s;*-/ (  S  AX-i.  )  ) 

RETURN 

END 


r.tNC"  ION  VMC  (  /  -2,  J  J  >1, 3J12,  NT) 

CO«i:»  < / 7 x  vi^/Glj  i/  N  j  ? 

V  No'Pr  =  G1  *  ■*  i  A  N  (S3'7(  J  *?*  *  ?  -1 , )  /SD-aTAkISD^I  (  V^2**  2-1 . )  ) 
exn:  =  j:~-7n'j2 

r  ~  7  j  n 

E>0 
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TABLE  A-5 

SUBROUTINE  FOR  FACILITATING  ITERATION  PROCEDURES 

}  J  ° :  '  .1 1 1  s  .  y. r  n  l  i  i 

a  c  .j  <r t*  ,  o*  a*  j„) 

u  u,r,«L  :•  i  f  *  us 
r  <  --n  a  l  -  >jk  com 

iTfY  =  a 
lTr  —  =  1 

XT  =  XI 

5  n  t  Fuse  \<X1,S.33.  ,r:.Y;,ir. 

?  F  (  £  Ji(cI  )  ,  S  T  ,  0  C  N  V  •  3  )  GJ  r  )  :  n 

6  r.;OT  r  xr 

FI  =  FUM:,  M  (  XI ,  rLD-»:  ,[>[.!  S 

STn'jc  =  ,  ~  *  U  ” . 

GO  v^<  5Q 

1 0  I  F  (  0  L  T  «  .  o  ~  . 's  ^  ,  n )  G"!  ID  i  g 

JrL  =  HE. 1ft 
XIP„  =  XI ♦DEL 

It-  (  XIPD.  i  i  .  f  .  0*.  XI uT  .  ?J  X  I  r  J  =  <I-L'CL 
f:?L'  =  FJ  >*C i  n  < X I ,  Swra_,  G:.rr  ,  Ay^  <0) 

Sl'jV'  =  :rn-c:>  /  (xjj  >x[ ) 

IF  IF  (SLOP  F.  ••  0.  ft,  0)  G  0  Tl  ’6 
Xi 1  =  XT  - 

If03j(Xl-Yil)lt:,i.i.-^5)  it  ’  3  G 
lF(XIl,Gf  .5)  ro  T  o  17 
X  i  1  =  <  A  *  X I  >  /  2 . 0 

G.1  T  0  20 


17 

if(xii.l. . j>  gj  to 

25 

ATI  -  (3  + >1)72. (3 

20 

I  F  (  I  *  f '  *  ?  .  G  l  .  3  )  SO  i 

0  2  r 

25 

XI  S  XII 

I  '  E  i  =  I  ”  ”  S'  p  .  1 

IF  (I  *  c.F'Y.uF. .  N  ^  I  Y  •  0  ) 

so  r o 

26 

IFir.-Y.l.L.O)  GO  TO 

30 

X  =  £ 

u:L  -  (t>-?)/2  0. 

FS  =  FUNr<,N<X,5w^r,"L.l  A,l) 

SF  o  -  SIUL (1.  0, FSt 
XS  =  x 

00  27  I  =  ? , 21 

X  =  X  ♦  DLL 

F  =  F  jf^r  N  <x  jS^pr,  ~rL  f  *  ,  I ) 

SF  =  GIGV(i.o,F) 

IF  (A^S(Sr-SPS) .ST. 1. C)  SO  Tj  21 
XS  =  x 
F?  =  F 

27  SF  S  -  SF 
GO  TO  3C 

28  XT  =  (XS*X)/2. 
i^cFe  s  2 

IT  F  Y  =  1 


~  > NA  « 
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TABLE  A-5  (CONCLUDED) 

G ' '  T  "i  5 

30  x  =  a 

OrL  -  O-ij/20.P 

P'-TM’  35,  f  N P  L 

35  For (lHb.MP,  15-tJir  sjt  :j'Jv,e%G) 

36  ftHlTFR3  =  ,I?,5X,5-iXi  =  ,=Il3.bf5<, 

6HXI1  =  ,  E.  15. 5/1.1  -  1 3.6, 5X, 

3  ‘♦HA  =  ,wl3«b,5*,  4HB  =  ,I13.b) 

Ti  45  1=1,21 

r  =  rjN^:  \ (  X  ,f „OP_ , r: L f 4 f r ) 

31  I*-’’  40,i,X,F 

40  FO^^i’d-i  ,  5  X  ,  H I  =  ,12,5X,4HX  =  ,  _15.b,5X,  4  HF  =  ,bl5.b) 

45X=X*0rL 

^  DU”  =  C.  0 

it  A  3  J  j  =  •  F  £  L  c  b  • 

50  ’“T  U'N 

t 
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TABLE  A-6 

PROGRAM  FOR  CONE  FLO'J  CALCULATIONS 
(IDEAL  GAS) 

PROGRAM  rcNr  (  TnPUT  ,  OLTP  J  T ,  T  5  ~  :  F  =  T  *  °UT  ,  T  AF'  F  =  HUT  ='JT ) 

0 IHENSI OF  TMFTft  {?<  ^  ,  VO  (?  r  ,  nc  (  2f  3 )  ,  *C  ( 0  >  ,  PHI  (  2P  0)  ,  VH2  (2fl«>, WN? 
sr-t 7705  1 

1  C[“iC  (e>»*)OfLT»US'*jGA*'*  iT  u- T  A 
ir(rr>r(c||  q,ir 

1'  Gc=  ( G  l'u  -<  .)/?. 

w^TT"  (F,?ir,s^nr^)tjSC 

■>  (lhj  ,FFX,*TP\|e  FL"'*  OUI~ULA7ION*  ,//, FIX,*  (GAMMA  =*,F6. ?,*)*, 

1 ///,rlX,  »CPNC  HALF  iKGL"  =*,Ff  . *jFG~FFS*  ,//,5RX  ,*USUPF/C  =  *,Ff.A 
?,///,ulX,«rHFT3.nrq«tily).u/^.  ,cX,*U/0*-,  FX,  *  X/F* ,  4  Y,*FHI-QFG*,5X,*^ 
‘T2*  »/> 

th-ti  ( r )  =  r  f  1.  T 
7CT">  =  C. 

JC  (  n)  -UF'' 
nn  *  T=l,  iTO 

:kf  =  i  (nr  (i-t  n  *•? 

TFS  17  =  (  (  r-AM-l  .)  *UrF*tT(’-l  )  «•  (  CF*UCF*7S(T-1)  /  tan  (THFTA  (I  -1 )  /  FO  ) )  -G  A 
l*ntC  IT-ll’  (VC  <1-1  >  **?)-GF*  (7F  (I-l)**3>  /TAN  (TH  FT  A  (T-l)  /PCI  )  /  (  (  (GAM  •• 
21  .  )  *  (  7C(  T-l)  *•?)  /?  .  )  -FT*  ’JCF) 

T  Hf*  '  (T)  =  THFTA  (T-l  )  ♦  "1rH.;TA 
7C ( J >  -VC (T-l » ♦OTHFT A*CE  7 1 7/ 

ur  (T>  run  r T-i  >  f^thft a*  (T-d  /fdm cthet  a**?>  *oeriv/  ( 2.  *  (ro**?>  ) 

wr  (T)  =SOf  T  { (UC  (I)**  2)  ♦{7C  C)  **2)  ) 

7-?<I)=S  ",ft(  wr  (  j»  **?/  (G"»  ( 1  .-VC (I)  **2)  )  ) 

OHT(T)sA’AN(Vrm  /  UC  (I))*®')  4-thETA  (I) 

CTT-,-.f*{Tf*l  (°  H  T  (I)/'1!')  Ml  .  /T  A  N  (TH"TA  (I)  /PD)  M  /  (  S I N  ( 2 .  *THET  A  (I) 
1/RP)-TAN  (fHl(T)/R3)MGA  ‘♦HS(  ?.*TH"TA(  I)/RD)  )  )  )  *  (SIN  (THETA  (I)  /RP>  * 
2*q>) 

W*2  (T>=S°PT(  (  <«GA**1.»*»  2)  *  (WMSTT**2>  /  ((FIN(TH;TA(I)  /RC>  )  *»2)-«*.  * 
1  (  W'*1«-TT-1  .  )*  (Gf  m*WM  1  STT*1  ,  )  )  /  (  (?.*GAH*  HMlST  T-GAH+1.  )  *  (  (GAM-i  ,  )*  WHl 
2STT «■:>.))  j 

IF(W?(I) -V*2(T|)  4  ^  4  jr 

F  WRIT"  (S,F )THFTA(I)  ,7C(t)  ,UO(I),WC(I)rPHI(T)  ,7H2(I) 

F  F0PMAT(1H  ,FL7.2,FlL.t,?F«.L,Fq.?,Fn.T) 

■>  C  0NT  T  NUr 

t  THFTAS=THrTA  (T-l)  OTHETA*  (XI?  (I-D-VM?  (T-l)  )  /( (WM2(  1-1)  -7^2  ( I-D)  - 
l  (WM',(l)-vr?(i) ) ) 

7H?S  =  WM7(T-1>  +(IHFTAS-TH“TA  (1-1)  >  »  (WM?  (I)-WM2  (I-D)  /  (  THFT  A  (I)  -THET 
1A (T-l)) 

F*HTS=CHI  C-l)  ♦(THe-'rAS-T  (£TA  (T-l))  *  (Pril.(I)-PHI  (1-1 ) )  /  (  THET  A  (T)  -THET 
1  A  ( I  - 1  >  > 

s A= (PTM ( (THET AS-OHIP) /R 1) )*«? 

S"= (rIN(7  HFTAS/PO) ) **? 

sr  =  v'^s**  ? 

7H1=F'T'»T(<9.*(GAM-1.)*S0*SA)  /(SBM?.*OAM*SC*SA-GAK*l.)>  ) 

MRITr  I6,71VM1,THFTAC 

7  FORH»T<lHt,l*qxf  »M1  =*,F’,L,5Xf*TH€TA  S  =♦,  F7 , 3,  *OEG*) 
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TABLE  A-6  (CONCLUDED) 


7  tusrm*?/<cfrmi\  ,-cjsc»*2)  j  n 

TT.SH<=PCSf*f*  ( (G4*-  1, )  *(  -.'"l**?)  *5?3*2.  »✓  (  (GAM*l . )  *  ( \Hi**2)*Se> 

Tf'C'=  (  Ga-/  (GA*'-1.)  )  *  flLCC-  (T^SrlO  »ALC G  ( t*~  SHK) 
rF  =  (S»G'-»sr)/(i1  »rF*{  7m?3p«»?»  > 
oc'ccc-fr^wx*  <&6>V  {G6M-1  .)  )  * 

T  £^S>f-TP<~  (-•«■*  ywF 

U*IT"  (f,  «>rRS«<,  rr,<:PF,  ”SHK,~PSPF,DS' 

FflpM*-(l^<-,i,iyt»(F?/Fl)  =*  ,r?  .  (^a/Pl)  CURF  =*,F7.3»/> 

1  **»*,*  (T?/TDSH<  =*,c7.?,13y,*  <T2/T1)  St'pc  =  *,  F7. 3,  //  ,f  lX,*ns/P  =  *,F 
2a.^  ) 

GO  1 
STOP 
c  rgr 
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TABLE  A- 7 

SUBROUTINE  FOR  CONICAL  SHOCK  CALCULATIONS 
(IDEAL  GAS) 

jl  •-  /w:  ?  ,0  Lil.f.  ms,  i  -it  S  » T  Hr  .  a  ji,  JS-.S,  VH2CS,  PHICa,FRC 

1 .»  ,T  A  OS#  7-* ;  S;> ,  P -sc  h  <  ,  ~  c  S  ii<j|  ’  -  S  S  i  V'  Hi?  z  kF  i  9  F  '  S  PF  S  *  T  S  1<FS ) 

r  yf  r  ^  *  s  i_  r  vi'  l  ~ 

“  V  1  ^-VELT,}^,  VM1,  TH£  :  ,  OS  \l  IcZ  ,PMl  -  ,  P\C,  TRC  ,  GF,^K  ,  7H2S  ,?  ' 

!  ji\,'  7S  HK , °HIS, VM2S-  F,»  -  S"-|*  'o'  r 
w  ‘1  =  7V1S 
j  :  l  7  :  L  T  c 
o,,,-i  =  "  SMS 
r  H-  =*  M£  C 

~  .29:  ■’791-i 

,  •‘■A- 7  1  *  ( 1  (DEL  7/  .)  -  (  (0  £u!  **£  t  /S  J  uT  .  )  ) 

j'3"  G  =  ,<j  (G^*  (  7  MG* *2)  /  (i  ,  +  &F*  (  VMG**2) )  ) 

'■i  -  (.  9Lc-.  r  '  9*DrL"  ) 

JCS-' N=S7,;  1  (GF  *^NS/(1.+  jF»J1Si)) 

US'  "V  -S  7t’  ~  IGF*  ( Vil**  2)  /  (1  .  *  (VMi**2>  )  ) 

:  *-L  L  *‘b  vi~  ~  (USC  ,  US.  G,  UiC  U,  Jj  ;-X  ,  .  '  Jolt 50,  .  „  ?i,  :GHCUCLC,LMClC  ,  X) 

I  HE  FS  ='.  HE  ras 
OS-  S-  OS"1 
4  m?  C  ~  =  7  *?r 
oh:  l  i  =ch'  ~ 
o : :  :•  =  f°  r 

7  *i  Sc  =  7  '•*  ?  _ 

3i’LF<c=  ’>  SHK 
SH< 

“*  n  :  S  A  rPHI'- 

FS  =  7F‘2$*F' 

-r  cr  r S==RSRF 
T  i'sS  F  r  S=  T  •  T  KF 

cr  UA  N 
ENO 


RJNfT  iOH  lNLLC  CJj''  ,  3UM1, 0UH2,  NT) 

COMMON  THf,PfLT  ,3rh  ,  VM,  THEtaS  ,  GS*,  V.12C  ,PHIC  ,  PRC,  TRC,  GF,RD,  VM2S,PR 

lSMKjTtSH'-'  ,HHIS,  7M2S*  :.F,P-iSriF,T.  SRF 
jirtt  JSirtN  THETA  (9  .  “>  ,  VC  (91  J)  ,UC  (9  J  3)  ,  WC  <9.l)  ,  PH  I  (9J  u)  ,  VM2  (91  J  ) ,  W12 
1  ( 9'  r  ) 

TH:T'  (r.)  =  OEL T 
VC  ( 11 )  =7  . 

IK  ( 0»  =USF 
OTHET/'s.l 
N=" 

00  1  1*1, 9QP 

ucf«i  ,-cur  (i-i) ) **2 
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TABLE  A- 7  (CONCLUDED) 

- 1<*  =  c  j  <  :-i)  ♦  t of*  j  i-u/Tmh (Tii rs  <i-d/k 

1  1 >  *  (VC  ( T  -  '  »  *  •  7 )  -  ,  - *  (  J  (*-i )  ♦*  3)  /  ~_  .('-’Lf.'.L-:)  /Ki)  )  /  ( 

71  . >  *  (  VC  (  J  -1  )  *  *  7  )  /  -  ,  )  - C-r  '  JC  r*) 

•  v  *  (1  )  -  ,  r ,  £  "  j  (  I  -  ;  )  *  ’,Ti  .  Ti 

r  f  ( r.)  . ,  i  , c  « 

r  IFC  TM.- ■  -  :  (I  J  J  -  ,  ,  t 

•I  -  - 

v  '  ( I )  =  v c  r  -i )  +  21  i ~ :  *  c L  I  v  / 1  1 

j  ...  ( i )  -u:  +  j«-('rnETA'»*t)-‘D^Riv/(£.#(Fj**c 

(  :  )  -  S  2  -  ,  (  (Ur  (I)  •»*_  )  ♦(73(1)*'  2)  ) 

J  (i)  =  57'T<W(.(D**.-/(G-Ml.-i'C(I>**2))) 

1F<T-D  ,%~ 
fc  T  f-  (  v  •  ?  (  I )  -  7  *  1  )  7,  5 

W'T  T  (b  ,  4  )  V^i  ,  \J*L  ( I ) 

f o- nr. *  ( 1-  ,irx,*wr-  *i  =*,-*.:*,*  isj^F  =*,"6.3,*  is  imhjss: 

Gl  t->  2 

3  Hi  (v  )  =  *’  AW  (  Vf'  (  I)  /IJ  ;  ( I)  )  *kD  *  HLt  A  ( I) 

A  ."11  S"  7  =  (?  .«•  (T  «N  (PHI  (  i>/  .3)  (1  ./TAN  (THET*  (I)  /PO)  )  )  /  (STN(2.*TH 

1/0)  -  H  ( I- Ml  {  I )  /-<})  *  (bA  «•(.  3  J  it .  *Trt£T  A  (I)  /  TJ)  m  ♦  (  SIN  (  TMtTA  (i, 

2  *  <  i 

W'9  r  )  =  S7  T  (  (  (  ( G 1 v  <•  i  .)*'  ?)*  (  *  r.  1  h  T 1  *  ♦  2  )  /  ((CIN(7HET4(I)  /  R  C )  )  *  * 

1  (  W-  li  TT  -1  ,)  ♦  W’liT'  1.  )  )  /  (  (2  .  ‘  bAH*WH13TT  -3  AH*  1.  )  ♦  l  (bAM-i 

23  ’*  +  7  .)  )  ) 

IF(  w-F<  I)  -')►  2(D)  u  ,  c 

~  got  o : 

r  H"  S  =  T  hlT  A  ( I  - 1 )  «-Q*  HtTo  (  H4?  (1-1)  -V12(  1-1)  )  /  (  <  WH2  ( I-i)  -VM2  ( 
1  <  w •  ■  2  (  I )  -Vi  2(D)) 

Vf  2  b.*  w-  2(i-l)+(HtTHo-THt’A(I-l))*(WM2(I)-WH2(I-i))/(THETA(I 
1*  (r-‘  )) 

°MiS=-HI  (I-i)  *(tHET.4  o-T  (ETA  (  T  - 1 )  )  *  (PHI  (  I)  -  PHI  (1-1)  )  /  (THEiA(I 
1  A  (T-l  )) 

3  A"  ( '.■  JN  (  ( '  H  E  T  A i - 3 H  I S  )  /O )  )**2 
L'  '7=  (c  IM  (  i  HFT A  h/\j)  )  *  *2 
3C=  VH  ??**? 

CMoU'-SZtT  ((2.MGL‘'-i.>*SC*3A>/tSoM2.*GAH*3C*SA-GA'‘Hl.n  >-V 
2,  (  VM1*-  2DR9-.D  *  +  1.  )  /  (G4M  +  1.  ) 

T  cbh<=P  =  S  HK«  ((bAfl-l.)*(V11**?>*S8*2.)/(<G4M*-l.)*{VHl**2)*S>3) 
OSRMGA  i/  (GO-1.)  )  Cb(TkSHK)  -BLOG(PRSHK) 

VV2C-  VP  2 < N-l )  ♦»  VH?  (.V  )  -V  *2  (  V-l  )  )  *  (T  H£-TH£TA  (N-l)  )/  (TrtFTA  (N)  -T 
1-1)  ) 

P  Hi  C-:  =H  I  (  N-l )  ♦  (PHI  (N)-FHT  (N-l))*  (T  HF-THETA  (-J-1)  )/(THETA(N)-i 
1-1)  ) 

V«- =  (1.  +GF*SC)  /(l.  ♦■u'r*Vl2C  **2) 

PkC=cRSHK*  fVMF**(G6  7IG51-1.)  )  ) 

TPCsT  RSHK*  V^F 

7*23' F=S7kT (U^C^*2/ ( GF* ( 1. - JSC**2) ) ) 

VnFS=  <1 .  «-GF*SC>  /  (i  .  ♦GF*  VH2SRF*  *2) 
pRSPc=PRSHK-»  (  VHFS**  (bAP/  (GM-1.)  )  > 

I  K5KF=TRFH<*VHFS 
2  RFIURN 
F  N  j 
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air  FORCE  AERO  PROPULSION  LAB  WRIBHT-PATTERSON  AFB  ON  F/B  *0/* 

CALCULATION  TECHNIQUES  FOR  INVISCIO  TWO-DIMENSIONAL  SUPERSONIC  — CTC(U) 
SEP  TO  MB  BCR  AS  TEN 

AFAPL-TR— 70-1023  NL 


AFAPL-TR- 79-2023 


TABLE  A-8 

PROGRAM  FOR  ILLUSTRATING  THE  USE  OF  SUBROUTINE  CONSK 


program  s  (in^jt  ,  clt-dt  , rA  -  fp  =  i  xp>jt» ta-5e  i soutp  jd 

2  ALL  (VM1,  J£LT,  G^,TH“,TmET  AbjDik,  v/M2l  ,  PHlof^Xu,  TK«“  ,  V^Sf^Srt 

1<)  V  M  ?  5  K  F  *  P  F  S  •  FjT-iSfFf 

w  vT'  -  (b,  2IGm:',,  \M1  ,  DrLT  ,  THET-  i ,  DSn.,  THE,  rfM2C,  ?HIC ,  PRC  ,  TRG,  VM2SRF  ,  P R 

1  S  ►  F  ,  T  F  S  P  F  ,  V  *■'  2  S  »  PH  I  S  ,  r  h  S  "t  K  «  T  <  >  r  K 

FrwH-,7(lrli,  =  4X,»CO.\<r  F^TW  C  AL  ■'  UL  A  T  10^*  ,  /  ,-j  8  <  ,  *  (GAMMA  =*  ,  Fo.  3  ,  *  )  *  ,  / 
1/,l7x,*mi  =*,FF.T,l?X,*TrLr a  --,Fu.2,*0£G*,//,5bX,*THETA  S  =*,F7.3 

2  ,  ♦-'EG*,// ,E*X,»L  3/x  -*,ci.4,///,55X,*»*T  ThETA  =•»  ,  F7 . 3 ,  *  DEG  *  *  , // ,  6 1 
3X,*W2  =*,F6.  ?,/,3?x,*r4r  -*,f  •  .3,*3EG*,/,‘iA<,*P?/Si  =*  ,  F7 . 4  ,  /  ,5  5  X  , 
'*  *  T  ?  /  ’  1  =*  ,F7  .*.,///,  f  •  X,*A’  3’J-F^Cl  «*,//, 61X,*m2  =* ,  Ffc . 3 , /  ,  ?  ciX  ,•  F  It 
iPl  =*  ,F  7.*.,/,C3X,»T£/’-i  =*  ,  F~  .,,///,  olX,MT  SHOCK  **,//,  bl  X  ,♦  M  2  =*, 
oFc.  T, /,?  c  X ,«  FHI  =*,F-.?,*- £,J  ,/,5iX,‘P:/Pi.  =*,F%*,/,58X,*T2/T1  =* 
’  ,  F  7  .  -  ) 

GO  TO  l 
END 
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TABLE  A-9 

SUBROUTINE  FOR  THERMODYNAMIC  PROPERTIES  OF  EQUILIBRIUM  AIR 

SUBROUTINE  THAT'  (  r  ,  F  i  ,  ,  S  *  ,  F>HO  ,  Art ,  C  3  ,  o  V,  -i  ,  A,  <1 ,  K£J 

T  =  Tr 

r>  -  n  » 

H  =  H  " 

^  -  c 

IF  (  <  «?-  1 )  1,1,2 

?  re  =  .f  F5r*GA*  / <  S  '.  *-l  .) 

IF  (<1-1)  3,4,=; 

7  n  =  c  r*  r 

3  =  I  .6  3  +  Cc*9L0G{T/5Cl'»)  - .  )ojr  *  ^Lt1  j  (°) 

GO  F 

4  '  r  'VC'' 

0  =  1  ,  e.  ■* *  r  o»  r l  o 0  (  T  /  L  0  0  .  )  - .  )53r*3LOG(a) 

or  "o  o 

f  IF  (<1-3»  ■’,8,a 

7  I  -  *  r'f ,« -:xx‘  ( (S-i. «.?♦.*  3#;»*3log(f  ))/cp) 

H  =  c °*  ~ 

GO  T0  f 

4  ►  =  "XX  (  (1.  F  FtD,  (7/5' ..  ,)) /.efts';) 

-  -  r*  r 

C  0  *  0  F 

9  T  =  H/CC 

p  =  r  XX  (  (1.63-i«-?P*  BLDG  (T/5  :  .))/.0665> 
f'  r  J  =  23. ar 

C  =  C®  -  .  "  bF  5 

F  HO  =  P*2l15.*A«/(15'~t>.*T) 

.1  =  S7RT  l32.1/4*&A“'*l3-5.‘r/.,-M) 

T  r  =  T 

?  f  -  3 

HC  =  H 
S r  =  ? 

GO  to  u?_ 

i  if  <><i-i)  ic,  11,12 

lr  CALL  FO  MS  (T,P,H,S,  =  HC,1M) 

H  t?  =  H 
it  =  c 

r-n  T0  13 
11  Tift  =  H/.24 
T2 a  =  H/.7 
90  14  i  =  1,50 

CALL  fqNc  (T1A,P,H1A,S*  RH  0 ,  A  M  ) 

CALL  EQNS  (T2A,p,H7A ,F,?hO, AM) 

IF  (A9S  <  <H-H1A)/H)-,G0  T0U)  15,1?, 15 
l”  IF  ( I  -5  ^ )  17,13,17 

ie  WRITE  (5,19)  TlA,T?A,HlA,H?A,H,f*,Kl 

19  FORMAT  (1M-,  22HTHA  IF"  Dn  NOT  CON  YE  RGE  ,  1 P5E16 . 4 ,5X  ,*HKi  *,I2) 
GO  TO  4 2 


i 
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TABLE  A-9  (CONTINUED) 


1'  TlfiN  = 

T2A  =  T2-' ♦(  (’i  <*-T2A)  *  (h-HZU  >  /  (SH.-H2A) 

1  *  1 =  *  1 '  r 

1  r  T  =  -  ?A 

i  .  =  c 

on  t  i  i? 

12  :c  i,?’ 

?r  T9  = 

T*>  =  ~rrr  . 
jn  " ;  J  =  1 , '  G 

OIL  l  F'MS  (TA,P,H,5*,*H'‘,A  9) 

1.1LL  rONc  (Tc3,D,H,  ,PH9,  £.•*) 

IF  (^s  (b-SA) -,n  ,rr'l)  2h,2F,2S 
?-  IF  (  J-FO  2t  ,2?,?7 
p  -  W  ■'  ’"’7  ( P, 19)  9, 59  ,  S  3,Tfl,  tp,  <i 

r,  o  T  •>  u? 

?<■  TAM  -  T9 

T.-5  =  TJK  (TA-T9)*  (S-S°)> /(31-FS) 

Ti  =  TCM 
?  u  f  :  T) 

T  ,  =  T 
Hi  =  w 
GC  T9  13 

21  ?A  =  ,"'ll 
D3  =  1. 

1C  ca  i-  =  1,5.1 

dLL  FQNS  (T,F4,rl,SA,PHO,M) 

GUI  FINS  (T,FB,H,S9  ,«ho,ah> 

IF  <A°S  (F-''A)-.QrCCal>  29, 2 9, 3 r 
3  IF  (^-5°)  31, *2, 32 

32  9‘?IT-  (6,19)  S,SA,P8,7A,»A,=3,K1 

GO  TO  1*2 

31  °  AM  -  PI 

•=3  =  rx  X  <3lOG(P3)M  (8L9G(:>9)-9LOG(P‘3)  )*<S-S9I> /(SA-SB)  ) 
? 0  P A  =  PAM 

2  4  P  2 

D  ,1  =  P 

Hr  =  h 
GO  -r  13 

22  T  A  =  H/  ,  2 A 

90  3T  <  =  1,50 
ofl  =  ,W>1 
FO  =  1. 

90  7i*  L  =  1,50 

CALL  LON?  (TA, PC  ,  H  A  ,  S  A  ,  ■■’HO ,  A* ) 

CALL  EQNP  (TA,PB,9B,^e,3H0,4M) 

IF  (APS  (S-SA) -.00  nnri)  35,35,36 

3*  IP  (L-5T»  30,3 a, 33 

33  write  (5,19)  S,SA»S3,TA,PA,°9,K1 
GO  TO  A p 

?7  °  AN  =  PB 

P9  =  FXX  (OLOG(PB) M (BLOG(FA) -SLOG (P3> ) * (S-S8) ) /(SA-SB) > 
3*  Pk  =  PAM 


90 


AFAPL-TR-79-2023 


TAELE  A-9  (CONTINUED) 

o  - 

cat  i  '  c  < :  u ,  ~ ,  ,  r ,  vh  j  ,  a  .*> 

Tt  -  ‘-/UC 


T  F 

(  '• 

jS  («TA-TL)/TA)-,f'O0C)l)  3  5,29,40 

L 

IF 

(x 

-?y>  ■,7,‘-i,4i 

-  1 

T 

1  *■ 

( S ,  1 

GO 

—  r. 

4  2 

7  *? 

T  f 

- 

TA-('A-t ")/?. 

1?  IF  (-2)  Lr,LZ,U; 

U  2  t  j  -  ,ooc<T 

T  ?  =  1.  np>: 

CfiwL  (t  1  ,?,-)!  ,SX,RHOX,  A*X) 

cau.  E  9  N  F  (T?,=>,H2,SX,x“l0X,a«A) 

C<~  -  (H2-H1)  /  (TE-T  1  ) 

:  'J  -  CP-i .9P6/A* 
bM  -  CF/fy/ 

t  =  '  7° T  (32.174*G"  ',*i5'*5.*T/a^) 
4  2  *"7 

E 


SU3ROMTINP  EOnS  <T,o,H,",ftHO*aM) 

T  =  T/l.® 

w  1  =  l.-FXX  (-227r../T) 

7c:  =  l.-EXV  (-339''  .  /  ”  I 

91  =  ?.*?.*tXX  (-113^0. /D+tXX  ( - io9 9 C . /T  > 
r  =  r'.f7.*EXX  (-2  2®  •  /T)  4-EXX  (-326./T) 45.*EXX  (-229JC./T) 

94  =  4.+i'.*FXX  (-227C0./T) <-6.*EXX  ( -4l 5 J 0 . /T ) 
ica  =  u  •  ♦  ill.*  f  XX  (-3860'’. /T)  4f  ,*EXX  (-5  3200. /T) 

9C3  =  l.*3.*Exx  (-’;  .e/T)  +  3.*EXX  (-130.9/T) *5.*£XX  ( -  22 1 J  0  . /T  )  *E  XX 
1  +5.*EXX  (-S79C 3./T) 

TQ1  =  ■?.c'*BLOG(T)  «-.ll-8|_UG(71)  ♦SlOG(01> 

902  =  3 .5*9lOG(T> -.-2-3L0G{ 72) 

=  2.5*310G(T>  f  ,3*n.0Gn3) 

=  0-  =  2.r*9lOG(T)  f  .•’♦9L')Gn4) 

303  =  2,c  *910G(T)  *■ .  74*.2*3lOG  (Q5A)  +.  3*810G(  058) 

=>Oo  =  2  • r.  *  3L  OG  ( T )  -14.24 

AK  =  EXX  (-59Cr’1./T«-2.*30?-30i) 

«<  =  EXX  (-1132G3 .  /  T  ♦  ?.*804-’,Q2) 

G*  =  EXX  (-166633. /T«-BQi*3}6-.  2*903-. 0*3Q4) 

ffi  =  (-.3«-GZPT  (. 1  4  +  .  e*(  1 .  +-4.  *P/ AKM  >  /  (2.*  (1  ,+4.*P/  AK) ) 

IF  (EA-. 1 99°)  2,1,1 

1  FA  =  " 

2  n  :  (-.4*F7RT  (,i6  +  3.8-,*(l.*4.*P/3K>l)/(2.*(l.+4,*P/BlOI 

Tf  (E 9-. 7 999)  4,3,3 

3  E9  =  r. 

4  Eo  =  1.  /SZPT  (l.t-P/GX) 

IF  (EG-  •  9C99 )  6,5,5 

5  EG  =  r. 

6  Zl=  1  .*Efi4.E9*2.*EG 
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TABLE  A-9  (CONTINUED) 

t  <*  s  2  3  .  ~  ‘  /  7  1 

?  =  1645. /AH 

'M0  =  /t»»-:x/  ( -2 7  7  c ,  /T  5 

/?-  -  f^9C./T>*cXX  (  -3  S  9C  .  /"  ) 

Ir  =  (  ???S  1. /T)  *  TXX  (-Ii39'1./T)  4-(i89P1./T)  *  E  XX  (-199^0. /T) 

G7:  -  (^t.  ,/T)  »£XX  ( -228 ./’)+(  326. /?>  *EXX  (  -  326  . /T)  ♦  ( 1140  0  G  .  /  T>  *  E  X 

IX  <-’ZoC  ./-)♦<  i-65  •  '  . /T)  *£<X  <-48o34./T) 

Q-~  ^  (?7-*nn%/T)*crXX  { -^77  j  J  ,/r  )  +  (2<+90  ,/T  )  *£XX  (-4151)0. /T» 

G5*r  =  (  7  36rfir  ,/n ‘EXV  (-336.  r  ./T)  ♦(3492Uf,./T)*:xx  (-5820C./T) 

)c  5t  -  (7il.»/T)»LXX  (-79.6/T)M9-+4.5/T)*EXX  (-158. 3/T) ♦ ( 11 JOCO . /T 

1  >  *  r  XX  U?rc,  .,/T)*EXV  {  -4  70  3  T . /T  )  ♦  (  3  7950  „  .  /T  )  ‘EXX  (-£79 

21!  r,  /*  ) 

HI  =  3. 5* VlP/Vi*;!3/ 01 
H?  =  7.5*x2r/«72 

h7  =  ?.^  +  C37/03 

Hu  - 

=  2.3* .2*C5fi^/9=£+.fi*Uc 33/058 

HH  =  ?.5 

Cl  r  -7Q14.H1 

c  2  =  ^0  2  *■  H  o 

S7  =  903»H3 

5-  =  j9“*H4 

££  =  nQ9  4-H5 

So  =  RQ6  +  K6 

XI  =  (.2-efl)/Tl 

X2  s  C.8-EB)/Z1 

X 3  s  (2.*EA-.4*E5> /Z1 

X4  s  (2  •*  E9-1 ,6*E3>  / 7 1 

X 5  *  2,* EG/71 

X6  =  YC 

TF  (Yl-l.F-ol  7,7,8 
x 12  G  =  0 
GO  TO  9 

8  x  it  G  =  9LOG(Xl) 

9  IF  (X2-1.E-6)  1", 10,11 
1 G  Y2LG  =  3 

GO  TO  12 

11  Y2LG  =  81PG(X2) 

12  IF  (Y3-1.E-6)  13,13,14 

13  X3LG  =  0 
GO  T0  15 

14  X3LG  =  3L  GG  (  X  3 ) 

15  IF  (X4-1.E-6)  16,16,17 

16  X4|_G  =  U 
GO  TO  18 

17  XULG  =  3L0G(X4) 

18  IF  (Y  5-1 .  E-f )  19,19,20 

19  X^LG  =  0 
GO  TO  21 

20  X rL G  =  BL  OG ( X5) 

21  IF  (X6-1.E-6)  22,  22,  23 

22  X  6L  G  *  B 
GO  TO  tk 
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TABLE  A- 9  (CONCLUDED) 


23  X6LG  -  31  nG(XE> 

24  S  *  \*  ( X  *  (il-XlLG)  «-X  2*  (SP-XZLG)  *  X  3*  (S3-X3L>)  *X  4*  (S*+- 
lXciLG«-.LC.‘*)*X6*<S'-Xf'LG)-6LJG(P)>/77d. 

r 7  =  ”*1.9 

rt  zF*Tp*(xi*Hi*X2*H2«-X3MH3M29E<3'!./T>>+X4*<H4+{5t>oJC 
12190  Xo’HcW'79. 

T  =  -r*l.« 

OHO  -  P*211f ,/<p*T ) 

RTT  IJ  iN 
E  MO 


niuTT  ION  PLOG(XY» 

IF  (XY-l.t-e)  1,2*2 
1  x  Y=  l.E-fi 
?  3L  )G  =  AL  °G<  X  Y) 

RF”  URN 
E  M"i 


function  r7PT  ( x  y) 

IF  (XY-i.E-*)  1,2,2 
1  XY  =  1 .E-9 
9  iZFT  =  SOPTCXYI 
RETURN 
END 


FUNCTION  E  XX ( XY ) 
SXY=FTGN(1.,XY> 
IF(A3S(XY)-8A.I  2,1,1 

1  XY=«:XY*89. 

2  EXX  =  EXPIXY) 

RETURN 

ENO 


X*LG>  »XSMSS- 

,/T)  )  +  X5MH5* 
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TABLE  A- 10 

SUBROUTINE  FOR  NORMAL  SHOCK  CALCULATIONS 
(REAL  GAS) 


1 

t 

3 


6 

7 


4 

5 


SuaUUTiNi  N'MSK  (  V-.l  ,  T  1 ,  »1»  •  /  "12,  f  2,  P2A,  GAMI) 

IF  (Giirtl-.l)  1,2, 2 
<2  =  1 
GO  TO  3 

<2  -  t  - . 

C  ft  L  _  MSh  <1  l,Pl4,nl,Sl,^Hl,  a  <11,  CV,GAM,A1,J,K2) 

VI  = 

M2ft  =  <-U*o.  *V11**2/(V41»*2*5.) 

DO  4  I  =  1,50 

P21  =  (31  ■<*  211b.  KHl*  <  VI*  *  2)  f  1  1.-RM1/KH2  A)  /  32. 17h) /211o. 

Ht  ~  rn+t+t+-+  It  Ml.  -  <  *H1/<M24  ’  **  2)  /50&62,  7 

C-„_  iHUa  (T 2, 32A , H2, S2, RHJ2, AM2, 3P»tV, GAM, A2, 1, K2) 

IF  (AaSU  <H02/<M1)  •  (  RM2A/R11)  ■  -.  UCJl)  5,5,6 

IF  (1-20)  4,7,4 

W-  I T  e  (&,C)  xH24,  t,Ml,k402 

F0-<*14'  (IriO,  1U<»334  (\hj<  010  >,3  1  OJNVtR&c - KM2A  = ,  lPt  11 .  4,  3X  ,  ?rt 

RE f  J  !N 

R42A  =  <M2A*  (RM32- Rrl2A)  72. 

V  r  =  SORT  11.-50062. 7* (H2-M1) /( Vl**2) ) 

VM2  =  V1‘VR/A2 

RET JRN 

£***} -  -  - - 
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TABLE  A- 1 1 

SUBROUTINE  FOR  OBLIQUE  SHOCK  CALCULATIONS 
(REAL  GAS) 


V]  ,  jri'.  .  .  -  ,  V*  *-  ,  r  ~  t  ,  f  ,  Ji  j,  T  o 

f  <  '  -  ,<£„  b  <r 

-'icr  j  vi  1 , ’  1 ,  r  i ,  ,  12»  ~  2t?2  t  sjArl*  ~.mO  i  »  Ml ,  .  i-;t< 

=  .  oi?<*r'  3 7 

r  ? ».;  =-t  l* 

W  -  1  =  v'  ■**! 
r  1  -  '  “ 

-  1  -  "  ^ 

»  •  ■*  r  »  •  -  -  T  H  ■ 

':.l-  '  rt*.  .  '  ("  1  *  -  1  ,  1  !  ,  .  1,  -  PM,  '  i,i,-  1,  1,  1, 

Till  I'V  L  ,  V*1  ,  ’  1  'll  ,  ~J  ,  "V  1 ,  T  H  «“T  A  ,  V  r  2  ,  ’  ?,  ) 

.  M  "  L  T  -  -  .  L  1  *■  )  1  ,  1 ,  5 
*  •  *  ’  1 
r  ■  -  -1 
4  ")=/’*  1 
r»  .  ’  J 

U  j=  (  (  V  i  ?  *  s>  I  M  ((’  Hr  L  )  *  \  T  )  )  /  (  V  1 1  *  b  I f  <  (  '  M  _  .  *  )?))*-,'  T  i 

.. ;  *  j  i  1 1.  ♦ t  r. 

ion  •  *  1) 

V’  j=  i  Z 
I  i  ~  r  2 
fia  ~  r  ? 

"l  H  t  I  =  ”  H  £  ’  '> 
t\  r  :  J  r  N 
Esin 


PJN:  ion  S<c  (LPfOwni,!  U“!2,  IT) 

-•  .1  pi  vt*j  k;  LT  j  >.  0 »  V ''  i  >  2 1  f3  d  t  &  •A  h  1 1  K  h  j  L  i  n  1  ?  T  -ii_  7  t 

T  HLT  -i  =  fcT  t  r.  (  (  (l.-J  -  )  ♦*?-  +  .  ♦UK*  l  (TAN  l  j;  lT  .1  )  ) 

i  j«.m  «n< j rkr»'a»>>^: 

p?= ( 1. ♦ (1 ,-UF ) *  SMI* ( *  ( S IN (  i"  H  -  I  A*  ^C)  )**2)*»1 

H2  =  <?  «♦(!  .-Ur  ^  /  c.)*  ♦(  (Sl  -Jt  J  1  f : T A ♦  ^ P )  ) 

:ut  THSir(T2,?^|i?>c2,<Hj2,"  2,3K,cy,uAi2,  A2,l,l) 

*  nor.  =.,hj1/?ho2 

s<:  =  jh -rh:° 

*  j*-  '-t 
en: 


, .  v  : ) 

1 1 , ,  i.  _ 


)  i  /  ( : .  * 

'  Z) )  * ^  1 
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«  TABLE  A- 1 2 

SUBROUTINE  FOR  ISENTROPIC  TURN  CALCULATIONS 
(REAL  GAS) 


VJ?r'',,.nr»  -  *•  r,  yc  r  !  ,  "l,  ’■ ’-l ,  ?"  rL  ,  : :  rL  t  ,r>-  V  LI,  vJ  2,r  r  ?) 
'r~.  '’17453"! 
n  r  cl  ”  -  ”  "  r  l  /r  F  r  L  *’ 

TFf  &  ’  r r  r  1 )  1,1,? 

1  K?~.l 

f,  n  -  •>  * 

7 

3  V“tll- A*T*’  (1  . /V-l»  /  " 

'•>  -  r-  i], 

T  «  _  T  , 

p?  ~r  1 
7’  •  -  ./"l 

•in  u  j=i,moc 
Pr  r  F<-  f  i  PL  T 

irn-.r.L  .^rFL » r,o  - ->  5 

CM  ~  H,'  ’  ’  (ZL,r!tu’  t~t:  -,r-A,  *  ,  I  r,'. V'»  '»K2) 

npit  ~ XN  (  "FI’,  V'*!-,''  f  ,**  ’  ,  ,  </*L,  V  1°,'  3,  1 -) 

7  r  =  7‘  ° 

T  *  r  1 

4  9*.  2P- 

5  V  U2=AST‘  <  1.  /V  '  ">)  f  ’ 

T?;*  7 

P  ?  =  r-  ' 

V  ?r7  a 
■  r  “  u  *j 
-  vp 
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TABLE  A-l 3 

PROGRAM  FOR  CONICAL  FLOW  CALCULATIONS 
(REAL  GAS) 

0  \  0  j  <  AM  RCONE  ( INPUT , OUi  PUT  ,  7  APE5  =  INPUT  ,  T  A0  E  6  =  J  J  TP'JT  > 

Ui MtNiluN  TaE  7  A (  2  00 )  ,42(200)  ,U2(20J)  ,*2  1200  )  ,  PHI  (  200  )  , 4M2  <2  0  J)  ,  k.Hj 
14(203)  ,  RHU1W  (2OU),T2(2uJ),02A(2O3),RHO.?(2OO),rl2(2jO),Vl<2UG),PlAW( 
1  2  3  0  )  , T 1W (200  ) 

G  - =66060.  13a 
RT=:>7. 29577951 
T.jJ  =  50n72.‘» 

1  R_4u  (5  ,*)  ;3EL7  ,  OThETA  ,USRF,  T  SkF  ,P3RFA 
IF ( c  OF ( 5 ) )  9,10 

lo  -KIT  E  (o,2)  OFlT, J3RF,TSRF,PSRFA 

2  FORMAT (lHl,5GX,*REAc  GAo  CONE  FLOW  GALOUtAT 10  4* »///>5iX, ♦CONE  H  A  l  " 

1  ^NGLl  =*  ,  Fb  .  ?,*3cjf<t-ES*  ,//,  57X,*JSURF  =  *,  F  7 . 1  ,* F  0S*  ,//,-,>  X  ,♦  i  o UR- 

2  -*  ,F7  .  1 ,  *  R*  ,  7x,*PSU-\F  =  *, 1PE10. 3, *AT M*,///,25X,*T HE Th-u£G*,3X, *42 
3-- Pi*, 4X, *U2-F0S* ,f X , *w2-FPb*, 3X, *0hI-jEG» ,:>X,*.M2*,7X,*T2-R*,7X,»0 
42-6TM* ,/) 

TnETA(O) =PElT 
42(0) =0. 

J2  (  Q)  =USRF 
02  A ( 1) sPSRFA 
T2 ( 0 ) =TS RF 

CALL  THAIR(T2(0),P2A(0> , HSRF , S , RH  02 ( 0 ) ,A4,oP,u4,GAM,A,y,0i 
HT  =riS^F+  (UoRF*»2)  /TGJ 
0  <HO  =  0  . 

00  3  1=1,200 

Ji\I4=-2.*J2(I-l)-42lI-l)M(l./TAN(TrtErA<l-i)/RJ)iORHO/RHu2(I-l)> 
FhETA  ( I)  =  T  H  E  T  A  (1-1)  OTHETA 
42  (i.)  =  42(1-1)  +  jTH£TA#l)ERI  4/ R3 

U2 ( I  )  =  U  2  ( I  - 1 )  ♦uTHETm*42 ( I - 1 ) /RD+pERIV* ( ( J7 H El  A/ RD )  *♦ 2 ) / 2 . 

W 2(1)  =  SZRT  ((U2(I)  **2)  KV2(I)  **2)  ) 

92 (i ) =HT- (W2 (1) **2) /TGJ 

CAcw  I  HAIR(T2(I)  ,P2A(I)  ,H2(X)  ,S,RHU2  (I), AM, CP, C4, GAM, A2, a, 1) 

4-12(1)  =W2  ( I )  /  A  2 

P9IU)  =A7AN(42(1)/U2(I>  ) *RO +THETA ( I) 

0 1HU=( RH02 (I ) -RN02 (1-1) ) *R0/ OTHETA 
UW2  =  W2  (  I)  *SiN  (  (THETA  ( 1)  -PHI  (  I)  )/■<□) 

Ul=JW2*TAN(THt  TA ( I) /?U) /TAN ( (THtT A (I) -PHI (I) ) /RD) 
41(I)=Ul/SiN(THETA(l)/R0) 

RHO 1  W  ( I) =KH02 (I) *  UW  2/U1 

PlAW(I) =P2A (I)  +  (RH02 ( I) * (UH2**  2) -RH01W( I) * ( Ul**2) ) /GP 
H1W=hT  -0/1(1)  **2)  /T  3  J 

CALL  THAIR(T1h(i),P1AW(I) ,H1W,SW, RHOW (I) , AM, CP, 04, 3AM, A, 1,0) 
IFUH01W<1)-RH3H(I)  ) 3,4,4 

v  WRITE  (6,6) THETA U) , 42(1) ,U2  (I)  , W2  (I)  ,PHI(I) , 4M2 ( I)  , 7  2 ( I ) , F  2A ( I) 

=  FORMAT (1H  ,F31.2,Fll. 1, 2F 10 . 1 , F9 . 2 , F 1 0 . A , F 1 0 . I, IPEIh.a) 

3  CONTINUE 

4  THET  AS  =  T  HE  TA  (1-1M-0  THE!  A*  (RHOW(I-l)  -RHOlW(I-l) )  /  (  UhUW  ( 1-1 )  -RHO  1H  ( 
il-l)  )  -  (  RHOW  (  I )  -RH01W  (I)  )  ) 

THtF=( THtTAj- THETA ( i - 1) ) / (T HET A ( I ) -T HET A ( I- 1 ) ) 
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TABLE  A-1 3  (CONCLUDED) 


J iIS^PHItl-l) ♦ iPHili) -PhI(I-I) ) *  T  H  £  F 

T2  S=  T2  C I  —  1  >  M  r?U)-T2(I-l)  )  *THEF 

p^5A=~2m  (1-1)  +  ( ~  2  A  (  j.  J  -p2 A ( I -  1 )  )  *  TH£F 

T  1  =  T  lW(x~l)+ ( T 1W(I) -7 iw ( 1-1 ) ) *  T  HEF 

?lA=->UW(i-l>*(PlAHlI)-PiaW(I-l)  )  *Tl£F 

U,LL  THAI*  <T1  ,r 1A,31 ,J1 ,RHD1 , A M , 3° , C , oA M , A  1, 0,  1) 

|/H1=  VI  (  I>  /A1 
^  <  IT  E  (6, 7) vrl , ThET A  3 

F  9  H  A  T  (1HQ,h9X,*H1  -  *  ,  F  7  .  t ,  o  X  ,  *  T  H  i  T  A  S  =  *  ,  I-  7  .  3  ,  *  J  lu‘  ) 

p\.jHK  =  P2SA/p  1A 

T:’SH<-T2S/T1 

P*  S*F=PSRFa/P1A 

T  Lb \F=  T *KF / T  1 

D3P= (i-Sl) /. 

^  IT  E  <  £  >  9 ) FPSHK, PRi ^F, T  kjH<,  T.R;,RF,  DSR 

FORMAT  (  1H[j41X»*(P2/P1)SH<  =*,F7.3,l3X,*(P2/pl)i>U<F  =*  » F  7 . 3  >  /  » 

U; X,* (T2/rl) SH<  =  *,F~.3,13X,*<T<?/Tl)*U<F  = * , F 7 . 3 , / / , 6 1 X , * j Sx R  = * , 
23.4) 
j  TO  1 
i  DP 
•O 
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TABLE  A- 1 4 

SUBROUTINE  FOR  CONICAL  SHOCK  CALCULATIONS 
(REAL  GAS) 

'  :i  '".J'J'  Vi-  L  Lb,  U  A  >.  ,  J-.  T  S  ,  T  -i:  T  ,3 

1  ,  7  ,  t  ~  y  ■  1  1  ^  ^  J>  ~  ^  ’  bl1'  J|  r\  \  jf  *  HI  >J(  .I’i^  j  '  f  ■  _ 

b<’LsU-L  •  "►T  * 

7  *i  “0 1,  '(1!)  .  '  HF  ,  „,l  *,«/'•:  1,  7  •  i  »  7  L.  T  v 

1  T  .3  H<  ,  3  r>  J  ■  ,  -  c  r .  -  F  ,  T  .•  „  F  ,  3  -  1  <  ,  7  \  3  r ,  -■  -  >  ►  ,  ■,  1 .  r 
J  ■'  i  -  J*  1 
T  ’■  =  Tib 

^  7=  :.r  l  ” : 

*  >r_  -7 

J  '  j  -  JK  1 M  !  .  -  ( 

7  F  =  OA  L  ’  *  (  v  *1 1  *  "  2 )  /  -  7  '  "  (  (  7  *  *  2  )  -  1 .  ) 

3  ■;>'<F  =  t.  +  ,j2-*'3r- 
V,'  1  I  -  1 ,  •_  •? 

*  SF.F  =  l  JF 

j:>  2  J=i  ,  -  ? 

)b-  FG--A.*  ,T(*  .i^K*~A) 

J  V-  F-f,=  lA  '  . 

J-^FMX  =  L  j  l’S  (  :  -5RP*T  A*  1  7  .  ) 

TALL  N  i  W  T  ^  (  J  b  F  ,  J  3  y  F  7 ,  J  S  Rr  -1  ^  ,  1 J  7  .  F  1 X  ,  1  .  7  -  -j  ,  ,  o  ,  1  „  .  , 
w  NC?-  '  - .  >F  -T  '  ~  -.r 
I  F  (  J-  r:  1 )  ■ ,  6  ,  r 

H  -  I  "I  -  It,'  ) 

HFMl'tH  ,*TcH!>  i’r  mT.J’I  21'  ’'J*  C  J'-WL  <••:>') 

IFl  A3S  <3 ,v:? ) -l.t-oi  3,7,° 

7  3 t p -T '  ~r 

■v...  “4  Ci  =°-  SF-, F 

7 r=  (T  i-.Ff  -  i  . )  /  .  oc7 
I F (  I  —  :  j  >  , 

X  -  I  T  ^  l  ,  1  3  ) 

FikUT  (tM  ,  ♦'>,[  I  r  F~:  A  7  I  )  4  Jij  4u  T  CjW-J,-- 

ifi  «3S(rcn:i )  -1.L-: )  ,1 

3  \S  F.F  -  a"  S  F 
T  :iL  TA  SS=  T  Hr  T  Jr 
J  j  \  S=  LS  f\ 

J  2CS  =  7‘-,2' 
r'-il  CS-',3Tr 
3  CS=r<C 
T  „S=T-  r 
J  j2SS  =  '^h^  ■; 

3  bH<s  =  r  i-bH< 

7  -,3  HK  S  =  T  P'  -K 
3  asS=3Hi'; 

X  32  bKFS  =  \M2>  RF 
3  Xb  PF b=  F  PF 
T  ~S  F  F  S=  TPS^F 
PrTU?N 
F  NO 
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TABLE  A- 1 4  (CONTINUED) 

rj.4f  T  '  .  .  ?  (J3<- ,  ♦DJI.',  I'l 

:  •'*;  ‘4,J4  j!  41.  )  »  '“t*  '-m  '  1  »  '  ^  f  1  ‘  '  •  »  *  7  •  • 

:  7  H<  ,  -“T  .  ,  •  :'F  ,  *  '3  f  ,  •  '  .  3;  >■ 

i :  ■  L ,, ■- t  r  i  * h (  * r  -z  : : ) ,  j •>  (  2  j  i )  ,  1 2  ( "  i )  ,  , :  (  7 : :  i  t  - 4 1  ( ° ' : )  >  i ■* -  :  -  -)  » '*  • 
1 1  ~  -  -j )  ,  *  ? !  - "  :■ ) , "  r  r "  ' ) ,  *  -vi  'i,  :  v.,  t  ■  ,,*»  4 

l” ) ,  -tni'i  < "  - .  )  ♦  *.  izc-D  j:  5 ) 

*  —  „  r  ■»  a  •  j  . 

^  .  *  i  — 

'=  -  ’ . 7  ■}  ■”  i‘ 4 

T  \,J  -  j  ,  ?  "  ’  .  u 
1  ■_!  '  r  .  ’ 

;r.r:,(4it  ■ ' .  T 
i  ■>  (  j )  -  7  . 
l  -(  r )  =■  iv  >- 
■>.’;())  -n  '  -  F  *  M 

r .  (  ' )  -  t  ,3  ^  f  *  t  >' 

:  „t. l  ;-u  t  :  ( *  ?  n>  ,  ,  ‘C  >f,  ..h  j.’  (  j;  ^  ) 

■r  -  n  i  •  r  +  (  1  j  ’  F  ■»  r  )  /  *■  o  j 

l  f  j-  I.  . 

•j 

; •->  3  1=1, 

j  ■  j.  )2  ( j.  -  j  )-<»?(  i  -i  )  -m  <  i .  /  T ;  ( i  -i[T  3  ( l  -  i )  /  ;* )  )  *•  w-io/  .h  j_  1 1  -  i  ) ) 

*  1  ,  (I)  =  'urT!  (1-1)  *'-H‘3Ti 

i  r (  n>  i, r ,  i 

I  F  (  T-U'-T-^Tk  (  r>  )  2,1,1 
■1=1 

J  M  I)  -1?  (  ’  -1  )  ♦  :)T  HtT  t.*n-E  ;  17/ '<"1 

J’’<")=J2C'-1)  OTHCT  A*  V?  ( 1-1 )  /"'  j  +  <T\/*  (  (  1 3  H "  A  /  •  3  )  4  *  ?  )  /  0 . 

?  (  1 )  =  ~  7  :  (PJ?(i)**2)  +  (\/2<:)**2>) 

HM  1)  =MT-(u2(Ij  *  *  2)  /"SJ 

I  F  (  h?  (t  )  )  -  ,  7  ,  <5 

•I'M  J)  =1. 

/•LL  ’■4-  T  ("MI)  (I)  ,  H?  (  I )  ,  3  ,  \  H  j  2  ( .1 )  ,  A  1,  =  ,  :  V  ,  3  A  A2  ,  « ,  1) 

1  ( I )  =  ^  7  (  T  )  /  £  ? 

-m;.).-;),  an  n  /u  (i))*«f’H:n(i) 

)  -M  r  -  ( -o  t'  (  T  >  -  J  2  ( 1  -  1 )  )  *•’. ) /  "1  * n l T  *, 

J.I?  -.-!•>  (  T  )•*■•:  ?i  (  (T  lET  f  (l)-Pnl  (I)  )/\J) 

J  i=  yw?»  t  (t  t  ( I)  /-  /  T;; ( 1)  -jh  r  ( 1 1 )  /  <:  > 

J  1  (  I)  =Ui  /rTNH  HK  Ti(  4.  )  /r-~) ) 

?  uj  i  1 1 )  =•  •  -i  u  ?  i : )  *  j  *  2  /  "  i 

M-W(i)=3’i  (I)  +  (\HT2  (I)  *  (M  4  2*  ‘2  )  -  <Hri  l  W  (  I  )  *(U1*,'2)  )  / 

•(  1W  =  HT-  (  VI  (.!  M  »2>  /TG  ) 

I  F  (  H 1  W )  ) ,  J  ,  K 
H  1M  =  1  . 

:  ^  l  u  '  h  u  i ( r  i  w  ( i ) ,  - 1 1  w  ( i ) ,  h  i  w ,  3 )  <  i ) ,  2  hg  *  ( i )  ,  a  i j  ,  g  a  m  ,  ►  <  i ) ,  i  ,  i  > 

I  F  (r-.H  Cl W  ( I)  -c  hOw  ( l)  )  3,4,  4 

OCNTIMJl 

-  ^rTiS=TH=T5  (1-1)  +DTHfc'TA  *e?HOW(X-l )  -vMJ  1  i,  (  1  -1  > )  /  <  ( InO  W  ( I-  1  >  -  *H  J  1  W  ( 

II  -1  )  )  -r’HOw  (  T  >  -- .H01W  (  I )  )  ) 

rtiLf  =  (THrMS-T  HE  14(1-1)  )  /1H3T4 
JHi  3=f-Hl  ( I  -1 )  +  (fHI(1)  -PHl(I-l)  )  *T‘iLF 
T2S=T2(T-l)MT2(l)-T2(I-l))*TlrF 
P2!iu:t'2A  (l-l  >♦  (F2A  (I>  -P2«(I-1)  )  *TH£;F 

>/'tlF=y/l(I-l)/ft(I-l)t(Vl(I)/A(Ii-'/l(i-l)/'(T-l))*T  HLF 
>  3  =  \/*il-7MlF 

I  (C^in  ( I -1  )  ♦  (T1W(I)  -Tiw  (1-1)  )  •-HEF 
3  1C  =F  1®  W(T-i)<.(fl„W(T)-°lA^(I-l))»THtF 
T  'FFr7?  (  i)  /t  AC 


TOO 
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TABLE  A- 14  (CONCLUDED) 


( ~) /ri~ 

7*2  =  F/a  <C) 

'  r.;jHK=T?5/ra 

y>?c=V**2<T-l)MVrt2(l»-V^2(I-in',T-li;.r 
•1  WI=SM  l-l)  M  SW(  I  >- (I  -  1)  >  •T-iCF 
.JS\=(L->''<I>  /.  Pfcrv 

T  -iE"J=  (THr-Iif';  (N-l  ))/(7McTm(  J)-T'«rT4  U-l)  ) 
V  "2C=  V“2  C’-L  )  ♦  l  N)  -\1«  2<N-1>  )  *  MIN 
JHTC=r^I  (Ni-1  )  ♦  (PHI  <  M  -FHKM-l)  )  *TM£N 

T  ?C=T2(f;-iM-(T2(N)-T^(M-i))*TtirM 

^?LA=F-’*t  (N-l  )  +  (►  ,?A(fv)  -P?M  7-1)  )  *THPN 
:~r  =t?:/t  a 
5 "  L  =.J2Cfi/°l 
Rt  TIRU 
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TABLE  A- 1 5 

PROGRAM  FOR  ILLUSTRATING  THE  USE  OF  SUBROUTINE  RGCON 


P  .vJGPA  1  XbCONL  (1NPJT  ,uUT  PJT  ,  =  ^  D  =  J  JT-’JI ) 

1  ^rA  0  (5,*)  v/Ml,TA,Pl„,  oALT.TH" 

I F( EOF ( i )  )  3, A 

-  CAlL  hG  C  0  N  ( \/  H  1  ,  T  A  ,  =  1  A  ,  j  *  L~  ,  T  H  t_  ,  T  H  c  T  A  S  ,  J  3  r< ,  2  J,  P.vi,  T  -  jc 

l-’.SHK,  7PSHk,PNxS,  V  rt?S  -  F  ,  P  PS  <F  ,  T  \S  IF  ) 

W  -ITE  (F, ’’WMt, TA,  =  15, OALl, THETAS,  IS  20, F  US, 2^- 

1,  =*SHF,  THSRF,  V^2S, PHIS, PRSH<,7*SrtK 
3  FQRMA"(1.i1,"0X,*^PAl  GAS  C'V4£  Fl-JW  C  AlG  Jl-*  7  IO.«*  ,  /  /  ,  o  1  X  ,  •  U  =*,r  j, 

1,  /  ,  ASX,  *T  1  =  *  ,  Ft  .  1,  *  ,*  ,7  X,  *  =  1  =  *  ,  l  =  t  1  1. -r ,  *5  T-*  ,  /  ,  J  JX  ,  *  JiluT  A  = 

2o  .2  ,*UCG*  , // ,5oX,  *THFTA  j  =  *  ,  r‘  7  .  3  ,  *j£  G*  ,  // ,  ?  3  X,  *  J  o/ H  =  *, re. -,///, 
33  x  ,  *A  T  THETA  =  *  ,F7. 3  ,  *  3E  G?  ♦  , // ,  G1X  ,  *  H?  =♦ ,  =  b .  3  ,  /  ,  i  )  X  ,  *  Phi  ^,F7.*, 
t*  TE  b*  ,/ ,F  tX,  *  =2/Fl  =*  ,  F’  .U, /,G6X,  *  T2/T1  =*  ,  F7  .-,,///,  6  0  X  ,*u  i  Sj.’.FA. 
GE  :*  ,//,f;lX,*^2  =*  ,Fu.  3,/ ,35X,*P2/  =  1  =*  ,  F  ’ .  h  ,  /  ,  3  sX  ,  *  T  2/ T  1 
6//,  1 1  X ,  *  A  T  SH0CKt*,//,61X,*M?  =*,Fo.  J  ,/ ,  ^  1  X  ,  *.»Hl  =  *  ,  F  7  •3,*0rG*,  /  , 
7SX,*=2/=1  =* , F7. A,/ , 38X, *T2/*1  =*,??. 4) 

GO  TO  1 
3  STOP 
EMO 
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TABLE  A- 16 

IDEAL  GAS  PROGRAM  FOR  PROBLEM  1  -  SWEPT  WHIG 


Pc 0 'i'Jh  k7nG(7NPU3  .  0  uTplJT.r**-  r5sTN®0  1  •  i  •  ►  ■  •  i T  ) 

\  t  h  i.  *  0  ,  i. 

pn=. L174S*7 
®',  =  211G. 

1  PE  AE  f  5  ,?)  *■  0,1  0,p0,  Yl,  Y?,ST  ~,1,  STC-2,  II  G3,S  IG4  ,  lJ  1  -  j  a  ,  At  Fh; 

2  F o *  ■*  T(l?t  <=>.  ") 

IP<F'*p<G)  )  4,  3 

3  r.Ul  r  wrt-  (L  t  M-OA  ,  mo,  Al  ,  GT“  1  ,  MFC  ,  i  Lf  E,  SlGE  1) 
c;t  L  SWft-  (I/'*' 90  A,*-  J  ,  A  If -(A,  SJG  ?,  MFO  ,  f  Lf-f  ,  S I GE?) 

CALL  >  WEr  Ml  ,  fit  £- -I  A  ,  31',  ■*,  *EU  ,  A'.!-:  ,  STC-r  7  ) 

CAlL  SwFp  (L^s'M,  I"1 ,  ALpM%  SIr-  4 ,  HE  u  ,  A  l  pp,  S IG-  4) 

JFLT-'l  =  ~:Grl-f  LF<r 

JrL  T  2  =  CT  Gc2*6LpL 
OF  LT !■»=-'.  IGM-FISF  T 
JLLTE  4=-? IGt 2-rTGE  4 
X  1  =  Y  1  /T  ftk'  r;ir.l*pn| 

*2=Y?/TAW  <SIG?*PD> 

XT=(Yi»Y-»«.(x2-Xi»,r?^{l?IG4*-<'))/(T4N(SIG3*F  .)  ♦  »  *  N  ;^.G4«r  4)  , 

<4=xi  *X3- y’ 

X=  X 1*  X  3 

Y  3=X3*TAv  (PXG3“F0) 

Y4=  x I  a  *'  <CTG4*RD» 

C«LL  I'GSX  (OELTFl,4-OfTp,PO,l.  4,0,L(ELTn,THEl»^l»Tl*pl) 

CALL  CGSK  lr)*fL'’E?»*r0fTC  ,  PI  , 1 ,  4  ,  P  ,  j  =  u  T  m  ,  T  HE  2  ,  r*  ? ,  T  l ,  c  l ) 

CALL  FXMLLLU3, 11»T1,!'1»  1.4,  VHU3,  1  *,  T  3,  p  x> 

C«LL  FXNl  rr.L  T  r  4,i2,T2,C2,i.  4,  ,  'A  4  ,  T  4,  p  4) 

FrPC*  (P3*  V34P4*Y4-:>1*Y1-D?-*V?  ) 

N=PC*  (  P  ?*  X?4P4*X4-p  1*Y1-F3»Y3> 

DPAO  =  N*?IMMLPHA*RJ)«P*COS(Ai.phA*Pl>> 

L:FT  =  N*ms  (AlFHA*fO)  *F*SIN(AI  rHA**9> 

WplTF(G*F)^0,P0,T0,ALP-'A,Yl,Y2,SIGl,SIG2,ilG7tSiG4,LAH3jA,M^0,SlG:. 

U,STGF2,SIGF3,SIGE«»,nFLTEl>TrLTr2,JFLTF3,0ELTe4,XllX2,X3,^4,X,Yi,v 

?4,0FLTM,‘'l,F2»kl7»44,Tl,T?,T7,T4lPl»t'2»p3»p4,OPAG*LIF7 

<5  FORMAT  <  H  1, 19X,*PRrSSURE  niSTPInUTIuN  nvpp  SwFPT  WI NG  *  /  3  3  >  ,  *  Of  Imp 
UNITE  ?PXN*,///,19X,*M0  =  *,F4.  !,3X,*P0=»,1PEA.?,*ATM*,3X,*TQ  =  *,'}Ppf 
2.0,*P*,7/3lX,*ALpH3=*,r5,l,*LFG»»//27X,*Yl=»,F4.2,*FT-,3X,*v2=*,r4 
3. 2,*PT*7/llX,»FIGl=*,FA.l,*DPG*f?x,*si&2=*,P4.1,*0PG*,3X,*STG3=*,P 
4  4. I* Jt G*3X*SlG4  =  *f  F4. 1 , * DE S* ,  // 30 X , * L AM3DA  =  * , F5« 1 , * OF G* , // / , 3 3 X  ,  * * 
5E0=*,  PG.3,//TY,*SIGF1=*,F5.2,  *Ot‘G*  ,  3  X  ,  *S  I  Gt  2  =  *  ,  F  5 . 2  ,  *  OEG*  ,  2  X ,  *  S I &F 
b3  =  * ,F5.2,  *OEG*,3X,*P7GP4=»,FG .2, *uEG* ,//5X, *0FLTE1  =  *, F5. 2, *fEG* , ?x 
7»*l)FLIt2=*fFF.2|*0EG*,2Xf»TEL'rE3=*’»FS.2,*lJEG*»2<,*D£LT£4=,fFG,2»*''' 
4FG*7/9X,*Xl=*,F5.?,*pT*,7Xl*X?=*fP5.2,*FT*,bX,*X3  =  *,F5.2,*P-r*,7X,* 
9X4  =  *  ,F5. 2,*FT*,//2qx,*T0TAL  L  FNGTH=* , F5.  2  ,  *  FT*  ,  /  /  2  7X  ,  *  Y  3  =  *  ,  F  4.  2  ,  *  p 
AT*,5X,»Y4  =  *,F4.  2f*PT*,//73IX,*DELTiAX=*,P5.2,*JLG*,//9X,*F‘l=*,PS.  ? 
1,  PX,*M2=*,F  <5.?,8Y,*n3=*,p5.  ?  ,  9  X,  *  *  4=*r  5 . 2  , Y  ,  *T  1=*  ,  F  6.  1  ,  *R*  ,  7 
2X,*T2  =  *,F6.1,*R*,6<,*T3  =  *,p6.  1, *R*  ,  3 X, *T 4= * , F&. 1 , *R* // 5 X , ♦ pl=* f lp 
3E9.3f*ATM*,2X,*P2  =  *>F9.3,*ATM*,2x,*P3=*,E9. 3 , *A TM* , 2X , * P 4-*  .  E9 . 3 , * 
4AT*4*,///,  2?x,*0PAG  =  *  »OPF0.  1,*L9S/FT  OF  WlNG  S  PA  N*  ,  /  /  2  2X  ,  *  l  I  FT  =  *  f  F  P 
5.1»*L9S/PT  GF  WING  SPAN*) 

GO  to  i 

4  STOP 
FNP 
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TABLE  A- 1 7 

REAL  GAS  PROGRAM  FOR  PROBLEM  1  -  SWEPT  WING 


P<ny-  AM  *  I  N  C-  (  T  N  F  J 1  ,  JUT=  JT,  r\=  l  5  =  I U  T  ,Tm  =0'J  i  -’LT  > 

MAl  4,  LIFT 

•  r)  =  .  1  174933 
D~  =  2 1 16. 

1  f  A  j  <  6, 2)  “  T ,  T  0 ,  c  0 ,  *  1 ,  Y  2 ,  SI  Sly  S  j.  o  2 »  »Io3»S.I64  ,  i_~>  "1  •  J  4  ,  ALT  4  L 

2  Pr,-MT  (12E  6.  0) 

IF(£0F(5)  )  4,  3 

3  CALu  SWE^  <LAMaDA, MO, AL =HA, SIS 1,-EO, AtF£,  SIG  .1) 

C«LL  i  WE  r  (L«fcjA,k3,.,l::HA,3l32,r.rG,fLPE,SI6£2) 

TALL  F  WES  (LAM?Sm,  Mu  SIS  3»  “E  j,  AuFE,  SIG-:?) 

CALL  S  Wr“  (LAA  c?A,1‘3,  A  L  ~  H  A ,  3IG4,M£o  ,  ALFr,STG  14) 

OL L T  :  l  =  F'GEl-A»Pt 

r'cXT  =  2  =  ST  G  F  2  ♦  ®  L  F  i. 
r» L T *"  j=“?  7  G  f  1-S  I GE  7 
DrL  1  l  4=-i  a  Gt2  IGF  «* 

X 1  =  Y 1 / T A  N ( ^TG1*5Q) 

X2=r?/T Aw (iiG2*RO) 

X3=(Y1+Y?«-(X2-X1)**A<(SIG  4*  ))/<*.  ANISIAS*  ,  L )  <•  T  A  N  (S  IG  4*  i-  3)  ) 

X4=Xl+X3-X2 
X  =  X  1  +  X  3 

Y3=X2*TAN  (5iG3*c.C) 

Y»»=X4*TAk!  ( jIG4»RO> 

LALL  rG?<<')FLT:l,*»E0,T0,P0,THEi,‘U,Tl,Pl> 
call  f'GSx  C,ELTE2,m-0,1  0,  P0  ,THE2,r12  ,T2,P2) 

CALL  -  G  c  X  F  (Tl,3l»'il,TtLrE3»?r  ,  0 ,  V"  J 1 ,  VMU  3 ,  r  0 ,  *  j ,  -1  X) 

CALL  ►GLXP(T2,P2,'42,OtirE4,2C.  ,  0 ,  Y  J2  ,  VMU  4,T  4,  Fn,  *  4) 

F  =  PC*  (P3*  Y3  +  F4*  Y4-T,l*  Y1-t*2*Y2) 

N  =  PC  *  ( P  2*  X2  +  F  4*  X4  -  3 1*X1-°3*X3 ) 

OF  AG=  6  *SIn(ALFHA*--')-F*:3S(Al  PH  A*  '  3) 

LlFr  =  N*COS  (ALFHA*c-*>  *F*3H(4L  FHA*  <  0) 

MA-  IT  e:  (6,c)"n,  FO,  TO,  Alpha,  Y1,y2,  SI  61,  SI  G2,  GlJ3,JI64,uA  ii  3l A, w-  0 , S * 6W 
ll,SIGt2,SIGt3,SIG£4,r.'(LTFl,JfLlt2,DrLT£3,0cLTt.4,Xl,x2,x3,v4,x,Y3,Y 
24,r-LT6,.-.l,M2,'i3,«4,Tl,T2,TX,T4,Fl,P2,F3,P4,"KAG,LlFT 
5  FuF“AT(lHi,iqX,*PVE*$lRE  OISl  hlttUT aON  O/EF  jW£PT  WING*/ 3  Jx , *  JF  iwF 
II  HIT-  SF1M*,///,13X,  *  MO  =  *,-'<♦.  l,3X,*Fu=*,  lPEi.  2,*AT  M*,  3X,*"  0=*  ,G-F5 
2  •  0 ,  *  >  *,//31X,*ALPHA=*,P6»l,*SLG*,//27X,*Yl=*,F4.2,*FT*,3<,*Y2=*,r4 
^.2,*FT»//iix,*SlGl=*,r4.1,*JEG*,3x,*SlG2=*,F4.l,*JEc*,3X,t3lG3=*,c 
44.1*'tG*3X*SIG4=*,r4.1,*OrG*,//30X,*LMM3DM=*,F5.1,*CEG*,///,33x,*- 
9£3=*,Fb.  ?,//7X,*SIG£l=*,F5.2,  *3  EG*  ,  3  X  ,  *S  I  oE  2  =  *  ,  c  5 . 2 ,  *  o*  ,  2  X,  *  S I  Gt. 

6  3  =  *,F5.2,*L)FG*13X1*Lj.GE4  =  *,FE.2,*4EG*,//GX,*Dt.Lr-l  =  *,F5.2,*FE&*,2x 

7  ,  *  Of  L  TE2=  *  ,  F  5. 2,  *ar5»  ,  2X  ,  *  K_  ’£  3=*  ,  F6 .2,  *">EG*  ,  2X  ,  *  ITU  I  4=*  ,  F6, 2  ,  *  0 
•‘EG*//9X,*Xl  =  *,F5.2,*FT*f7X,*x2=*,r5.2,*FT*,6X,*X3  =  *,'r5.2,*Ft*,7v,» 
3X4  =  *,F5.2,*FT“,//2'sX,*T0TAL  tENGTn=*,F5.2,*FT«,//(?7x,*Y3  =  *,r4.Z,*r 
AT*  ,5X,*Y4=*,F4. 2,  *  rT  * ,///3JX,  " DELI  1 A  X  =  *  , F  9. 2  »  *  JE  G*  »//9X , 1  =  *  ,  F  6. 2 
1,  9X,*H2=* >F5.2»8X , *  w  *  =  * ,”9.2 ,  3X , *M 4=*F 3 .2 , // 9X , *T i  =  * , Fb . 1 ,* K*  ,  7 

2  X , *  T  2=* ,  F  6.  l,  *RVf6Xl*T_3_=*,F_6.  1  ,*R*  ,  a  X ,  *  T  4=  *  ,  F  6 . 1 ,  *,T* // 5  X ,  M=*,P 

3  E  ■s  .  3  ,  *  A  T  M  *  ,  2  X  ,  *  P  2  =  *  ,  £9.  3  ,  *  AT  M*  ,  2X  ,  ♦  P  3  =*  ,  t.9.  3  ,  *A  T  M*  ,  2X  ,  *T  4  =  *  ,  E9.  3  ,  * 
4AT-*,///,  22X,*')RAS  =  *,  0FF8.  1,"  L9S/FT  OF  WING  SFA  W*  , //22X  ,  ♦  L  I  F7=*  ,  F  -i 
3,1,’LBS/F*  OF  WING  SPAN*) 

GO  T"  1 

4  ST  OF 

E  N  D_  _ _ _ _ 
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TABLE  A- 18 

PROGRAM  FOR  SOLUTION  OF  PROBLEM  2-  2-D  INLET 


.  r .  ,  > 


*  i  /  i 


*  )  •  . .  •  )  j 


■  s 

///*':  *  •.  « 


•  -  f  *  *  /  /  . 

'•  !  *  ■  t  *  1  i  v  T 

v  «  * 


•  ’  )1  ; 

,*  «  / 


•  ■*  -  >  v<  •  -  •  * 

.  v  /•>%/✓// j 

*  ■  t  •  •  J  *  •  •  ’ 4  *  . 


•  v  .  ♦*  i>  /  •  1  6  -  *  ) 

i  "*  1  •  ^  I  4  ^  !  r‘  •”  ^  ,  1  V  *  v  ) 

•  *  •  *  i  t  *  »  ‘  rH  /  .  i  v  -v  *  ) 


^  ‘  •  1  •  ¥  :  ‘  >**//*) 

"  »  V,-  -  I  (.  '  </  ■'.  • 


»l.v  . 


) 


•T  *  ;*  •/  >  /•  '  (•. 

' '  •  *  •  :  l  '  i  »  r>  l  ,  i  y  t  ,  \  V  |  4  f  . ;  ) 

'  '  -  i  i  »  •  .  <  I  .  »  ^  '•!»».  L  «.  i  •  I  i  ,  i  ) 

1  r  (  «  <  *  J  L*  !  i  .  «  i  i  *  .  %  \  )  *  ) 

»  y  *  i  /  *  !“*  *  3  4  T  :  *  r‘  J  *  ) 

*  1  ’  i  ‘  1  *  i  •  -M  •  *  ,  *  4  -  1  7  4  >  1  »  .  r*  4  C  *  i  ^  4  4  1  ) 

IT-  ;  .  ; 

1  »  !  •  >  i  -1  •  i  .j  j  _>  j  i_  T  i  * ' « >  t  /  J 

•ii  <  «  •  *  i  i 

‘it  i»4*>  4  r*  t  T  • 

*  !  ■  i  »  :  >  '  1"  i»  f  1 

I  ■  1  ■  !  •  i  *  v  »'  ’  * 

f  •  i  ! 

1  '  i  i  f  •  f  i  >  *1  |  f  J  ^  f  ^  I  ^  I  j  4  1  ,  ^  i)  1  4 

t  *  •  ^  ;  •  /  'i  • 

.1  1  •  *  t M .  i  i  Hi  ;  , r  '  i  I  i  ^ 
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TABLE  A- 1 9 

PROGRAM  FOR  SOLUTION  OF  PROBLEM  3  -  CONICAL  INLET 


F.-.JC  •'  A  !>■  i-3D  (INPUT,  C  UTPUT  ,  T  AP  £3  =J  liF UT  ,T  APE8  =OUT°UT  ) 

'FAL  K,:-l  ,M2 

r  E  AC  (5  *  2)  *C  »  r  u  ,  Pu  »£_£*"  ,C WL  ANG 
FOR  MAT  (i-cir.j) 

OF <E OFT*)) 4,3  - 

THE  =CE  lTMCWLANG-lElT)/2. 

CALI  C0N3K(MG»OEuT»l<  4  , HE  ,rHL7  A  ^  .HSR  *  Ml  ,PHI  ,  Pf  1,1  R  1 ,  VMI  3,  PS  SHK  ,  IF 
1SHK,PHIS,VM2SRF, PRSRF , IRSRF) 

*  l  =  T'"l*-l  3 
r  l  =  r  ,  1*"  'J  u 

CA.L  NChKfWI,~I,Fl,M2,  “2,P2,1«  4} 

tTAF  =<  F2M  (1,  *.2*  H2**2I  >♦  +  3*!>>)  /(PoM  (1.  +  .2-*  <MC  *2) >  *' 3.5) ) 
WtTTE(d,5)mj*PG»  ,C,DUT,CW1  A*J& 

FU'  MAT  <1H1,17X,*  LACK  P-ESSU.E  AND  PRESLU  t  U  CO  Vt  P  Y*  , // ,  27  X  ,  '  F  jF  C 
1  CNICAL  INlOTr  ,///,15Xi’rM0  =  *,'F4»i,AX,  *PC=*,  1PE8  2 ,  *  A', 3*  ,  4X,  •  *.  0  =  ’  ,  :  P 
2F6,  1,-  R  -,//,?  ftX,  ’DElTA  =«F-,,1,*CcG*  ,///,2bX ,*COWL  A  NGL  £=•  ,  F  5 , 2  ,  DEG 
3*“,  A<F//,2FX,*0jFAt  3A~  FCTUVIOJK*,/! 

H  ~  i  T  E  ( c  ,t)>'l,Pl,  i  1 
K= :*E <b H2,P2 ,  2 

rC '  MAT  <  1H  ,  13X, »  M1  =  *,F.#  2,  *.X,  ♦P1  =  '‘,1P£9.3,  :  ATM*  X,*T1=*  JP- T.  l‘P'  > 
FH-HAT  (1H  ,13X,*h2=*,FE.  2 , 4X , *F2  =  * ,1 PE  9. 3, *ATM* , OX , *7 2=*3PF7 .1***  ) 
W*~  'TE<c,4)THETAS,P2,LTAc 

FORMA**  { I'M  ,/T?BX,*7HrTTVS=^r9,2r*CtG'*,//,17X,*'‘~QUI?XC  3ACK  PC“SSU 
lRt  =  *  ,  lc'£S«3,+ATM*t//,l9X,*PPcSSURF  RtC  OVER  Y  <  PT2/FTD  )  **  ,  uPF  ft  .  -,////  ) 
C A!  I  *j  HA*T  CTu,FOVHH,£<r,RHTiA«OVCP?CV,GAM,AO*i3,11 
VO=MC*AC 

hts-hgm  vo^zt/bcg^.  ~  - 

CA.L  THAIR<TT3,PIC,HT,S0,RHT,AMT,CP,CV,GAM,ATiJ,4,C> 

XALt  V  GC'JT*  FPO  .CELTiTHE,  TH0TA3,T75P,tir,PHI  ,FK1 ,  T^I  ,VM2S ,  PRSHK , 

IT  •ShK,PHlS»i/l*2i,RFtF:<i.RFf 'RS'JF) 

UwL  neh;k( mi  , pi ,-m? rrz » pz  ,  c  > 

CAlL  THATP<T2,P2,H2,S2,RH2,AM2,CP,CV,GAM,A?tU,0) 

W  iiUS9)  -  "  - 

FJ-LATUH  , 2t*X ,  *,SE A..  GAS  SOLUTION*,/) 

w-:7tl  t;  ,6i pi, pi, rr - - ~ -  ~  . * 

HfuTtlS.nf  2,  P2  ,  T  2 

CALL  TH A  rr  (TT2,^T2,  (T,  02  ,TrFTT  ,  AMT  ,OP,T V  , G  AM  ,  AT3,  4 , 0  ) 

FT  Ar  =PT2/PT0 

WR  IT  E  <  6  ,  e ) '  HE”  AO  ,  P2  ,  t T A* 

go  t:  i 

STOP  ‘  - - - — — -  .  ... 

em: 
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APPENDIX  B 

CHARTS  FOR  CALCULATION  OF  LOW  MACH  NUMBER  SUPERSONIC  AIRFLOWS 

In  this  appendix  charts  are  included  for  calculations  involving 
normal  shock  waves,  two-dimensional  oblique  shock  waves,  swept  wings, 
isentropic  expansions,  and  supersonic  flow  over  zero  angle  of  attack 
cones.  The  data  for  the  charts  has  been  calculated  by  use  of  the  programs 
and  subroutines  discussed  in  Section  I  of  the  text.  The  charts  are  to 
be  used  only  in  cases  for  which  the  ideal  gas  relations  are  valid. 

Copies  of  these  working  charts  plotted  on  graph  paper  are  available 
directly  from  the  author  at  the  following  address: 

AFAPL/ROT 

ATTN:  M.  B.  Bergsten 

WPAFB  OH  45433 
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M]  ~  MACH  NR  PRECEDING  SHOCK 


Figure  B-l .  Parameters  for  Normal  Shock  Calculations 
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M]  ~  MACH  NR.  PRECEDING  SHOCK 

Figure  B-2.  Maximum  Deflection  Angle  for  2-D  Wedges 
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DEFLECTION  ANGLE/SHOCK  ANGLE 


AFAPL-TR-79-2023 


1.0  1.2  1.4  1.6  1.8  2.0  2.2  2.4  2.6  2.8  3.0  3.2  3.4  3.6  3.8  4.0 

M]  ~  MACH  NR.  PRECEDING  SHOCK 

Figure  B-3.  Shock  Wave  Angle  for  2-D  Deflections 
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Figure  8-4.  Mach  Number  Ratio  for  2-D  Deflections 


STATIC  TEMPERATURE  RATIO  ACROSS  SHOCK 
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M-j  ~  MACH  NR.  PRECEDING  SHOCK 

Figure  B-5.  Static  Temperature  Ratio  for  2-D  Deflections 
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~  MACH  NR.  PRECEDING  SHOCK 

Figure  B-6.  Static  Pressure  Ratio  for  2-D  Deflections 
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M-j  ~  MACH  NR.  PRECEDING  SHOCK 

Figure  B-8.  Dimensionless  Entropy  Increase  for  Large  2-D  Deflections 
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TOTAL  PRESSURE  RATIO  ACROSS  SHOCK 


EQUIVALENT  MACH  NUMBER  RATIO 
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Figure  B-9. 

Mach  Number  Conversion 

for  Swept  Wings 
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Figure  B-10.  Angular  Conversions  for  Swept  Wings 
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TO  DETERMINE  CONDITIONS 
DOWNSTREAM  OF  EXPANSION, 


Figure  B-1I .  Prandtl- Meyer  Angle 
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MACH  NUMBER 

Figure  B-12.  Static  to  Total  Ratios 
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M]  -  MACH  NR.  PRECEDING  SHOCK 


Figure  B-14.  Mach  Number  Ratio  for  Cones 
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STATIC  PRESSURE  RATIO  ACROSS  SHOCK 


AFAPL-TR-79-2023 


M]  ~  MACH  NR.  PRECEDING  SHOCK 
Figure  B-15.  Static  Pressure  Ratio  for  Cones 


122 


T,/T„  -  STATIC  TEMPERATURE  RATIO  ACROSS  SHOCK 


AFAPL-TR-79-2023 


Mj  ~  MACH  NR.  PRECEDING  SHOCK 


Figure  B-16.  Static  Temperature  Ratio  for  Cones 
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MpMACH  NR.  PRECEDING  SHOCK 

Figure  B-17.  Dimensionless  Entropy  Increase  for  Cones 
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